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Abstract 


The renormalised quark mass in the Schrodinger functional is studied perturba- 
tively with a non-vanishing background held. 

The framework in which the calculations are done is the Schrodinger func¬ 
tional. Its dehnition and basic properties are reviewed and it is shown how to 
make the theory converge faster towards its continuum limit by 0(a) improve¬ 
ment. It is explained how the Schrodinger functional scheme avoids the implica¬ 
tions of treating a large energy range on a single lattice in order to determine the 
scale dependence of renormalised quantities. The description of the scale depen¬ 
dence by the step scaling function is introduced both for the renormalised coupling 
and the renormalised quark masses. The dehnition of the renormalised coupling 
in the Schrodinger functional is reviewed, and the concept of the renormalised 
mass being dehned by the axial current and density via the PCAC-relation is 
explained. The running of the renormalised mass described by its step scaling 
function is presented as a consequence of the fact that the renormalisation con¬ 
stant of the axial density is scale dependent. 

The central part of the thesis is the expansion of several correlation functions 
up to 1-loop order. The expansion coefhcients are used to compute the critical 
quark mass at which the renormalised mass vanishes, as well as the 1-loop coef- 
hcient of the renormalisation constant of the axial density. Using the result for 
this renormalisation constant, the 2-loop anomalous dimension is obtained by 
conversion from the MS-scheme. 

Another important application of perturbation theory carried out in this thesis 
is the determination of discretisation errors. The critical quark mass at 1-loop 
order is used to compute the deviation of the coupling’s step scaling function from 
its continuum limit at 2-loop order. Several lattice artefacts of the current quark 
mass, dehned by the PCAC relation with the unrenormalised axial current and 
density, are computed at 1-loop order. An essential property of the renormalised 
quark mass being computed in this thesis at 1-loop order is the deviation of its 
step scaling function from the continuum limit, which was so far only known for 
the zero background held case. 
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Zusammenfassung 


Diese Arbeit befasst sich mit storungstheoretischen Rechnungen zur renormierten 
Quarkmasse im Schrodinger-Funktional mit nicht verschwindendem Hintergrundfeld. 

Als Grundlage der Rechnungen werden das Schrodinger-Funktional und seine 
grundlegenden Eigenschaften erlautert. Auch die 0(a)-Verbesserung, die zu einem 
schnelleren Erreichen des Kontinuumslimes fiihren soli, wird in diesem Zusammenhang 
dargestellt. Des weiteren wird erklart, auf welche Weise das Schrodinger-Funktional 
dazu dient, das Skalenverhalten renormierter GroBen iiber einen groBen Energiebereich 
zu untersuchen. Das Skalenverhalten sowohl der renormierten Kopplung als auch 
der renormierten Quarkmassen wird in diesem Schema durch Step-Scaling-Funktionen 
beschrieben. Die Dehnition der renormierten Kopplung wird dargestellt, ebenso 
die Dehnition der renormierten Masse, die mit Hilfe der PGAG-Relation iiber den 
Axialvektorstrom und die Pseudoskalardichte erfolgt. Die Skalenabhangigkeit der 
renormierten Masse wird auf die Skalenabhangigkeit der Renormierungskonstanten 
der Pseudoskalardichte zuriickgefiihrt. 

Breiten Raum nimmt die Berechnung verschiedener Korrelationsfunktionen 
bis zur Ein-Loop-Ordnung in Storungstheorie ein. Mit Hilfe der so ermittelten 
Koefhzienten wird die kritische Quarkmasse, bei der die renormierte Masse verschwindet, 
in Ein-Loop-Naherung berechnet, ebenso der Ein-Loop-Koefhzent der Renormierungskonstanten 
der Pseudoskalardichte. Mit Hilfe dieses Koefhzienten wird aus der bekannten 
anomalen Dimension in Zwei-Loop-Ordnung im MS-Schema die anomale Dimension 
im Schrodinger-Funktional berechnet. 

Als weitere Anwendung der Storungstheorie werden verschiedene Diskretisierungsfehler 
bestimmt. Die kritische Quarkmasse in Ein-Loop-Ordnung geht in den Zwei- 
Loop-Koefhzienten des Diskretisierungfehlers der Step-Scaling-Funktion der renormierten 
Kopplung ein, der durch die Abweichung dieser Funktion von ihrem Kontinuumslimes 
dehniert ist. Verschiedene Diskretisierungsfehler der Strommasse, die durch die 
PGAG-Relation mit unrenormiertem Axialvektorstrom und Pseudoskalardichte 
dehniert ist, werden in Ein-Loop-Ordnung berechnet. Ein wichtiger Diskretisierungsfehler 
der renormierten Quarkmasse ist die Abweichung ihrer Step-Scaling-Funktion 
vom Kontinuumslimes. Dieser Fehler ist in Ein-Loop-Ordnung bislang nur mit 
verschwindendem Hintergrundfeld bekannt und wird in dieser Arbeit mit nicht 
verschwindendem Hintergrundfeld berechnet. 
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Chapter 1 
Introduction 


Since the idea of matter consisting of particles was invented by the Greek philoso¬ 
pher Democrit (460-371 B.C.), there has been considerable progress in particle 
physics. Especially the invention Erst of quantum mechanics and later of quantum 
held theories have improved our basic understanding of elementary particles and 
their interactions. At the same time, large accelerators have made experiments 
in the high energy regime possible, giving large amounts of results to compare to 
theoretical predictions. 

The modern theoretical framework of particle physics is the Standard Model, 
which is a gauge theory based on the invariance of the Lagrangian under the gauge 
group SU{3)c X SU{2)i x U{1)y- The SU{2) x U{1) symmetry associated with 
the weak isospin I and the weak hypercharge Y describes the unified electroweak 
interaction itt 0 i- In this context, particles acquire their masses via the Higgs 
mechanism [Q. 

The SU (3) is the gauge group of the strong interaction affecting the hadrons. 
The idea of hadrons consisting of several point-like constituents was raised by 
Bjorken’s analysis of the scaling properties of the structure functions in deep 
inelastic lepton-nucleon scattering |^. This Bjorken scaling could be explained by 
the assumption that the hadrons consist of point-like particles called partons 
1^. These partons could then be identified with the so called quarks proposed 
earlier by Gell-Mann and Zweig in order to explain why hadrons can be classihed 
into multiplets P, |^. As a theory for the strong interaction, Fritzsch, Gell- 
Mann, and Leutwyler proposed a non-Abelian gauge theory containing quarks 
obeying an SU{3) gauge symmetry, where the associated quantum number is 
called colour [pi| . 

This theory, called quantum chromodynamics (QGD) and today widely ac¬ 
cepted as a model for the strong interaction, does however pose some problems. 
First of all, the theory is mathematically demanding, making it difhcult to get 
predictions for experimental results. Second, the basic particles of the theory, the 
quarks, cannot be observed as free particles in nature, due to the phenomenon of 
conhnement, binding them together in colour neutral hadrons. 
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However, despite of these problems, QCD together with the Feynman path 
integral approach ||TT|, turned out to be successful in the high energy regime. 


This is due to a property of non-Abelian gauge theories called asymptotic free¬ 
dom []^ 1^, 1^. Asymptotic freedom means that the coupling becomes small 


at high energies, allowing us to treat the quarks as free particles and to expand 
expectation values in powers of the gauge coupling. The coefficients in the ex¬ 
pansion may be formally inhnite, but can be made hnite by regularisation and 
renormalisation. Regularisation means introducing an artihcial momentum cut¬ 
off, and renormalisation means removing the divergences before removing the 
cutoff again. This method gave, for example, good results for the violation of 
the Bjorken scaling |^, and for QCD corrections to cross sections in 

electron-positron annihilation |jT^, ^ 

At low energies, perturbation theory cannot be applied, since the gauge cou¬ 
pling is large at these scales. For this case, Wilson introduced the concept of 
lattice QCD |^. In this approach, space-time is treated as a four dimensional 
Euclidean lattice, the inverse lattice spacing a~^ serving as an ultraviolet cutoff. 
Removing the cutoff then corresponds to taking the continuum limit. On the 
lattice, QCD can be treated like a classical statistical system, where expectation 
values of observables can be obtained by Monte Carlo simulations. However, 
due to the limitations of computer power, the simulations are restricted to small 
lattices, making safe extrapolations to the continuum limit difficult. In order to 
improve the situation, Symanzik has introduced the concept of systematically 
removing terms of order a both from the action and from the observables, thus 
making the theory converge at a rate proportional to rather than a P5| , ^ . 

In any case, both at high or low energies, the free parameters of the theory like 
the coupling and the quark masses have to be hxed by a set of observables. Since 
these observables have to be calculated in the theory, one has to take observables 
computable in the respective framework, i.e. lattice QCD or perturbation theory. 
This means that, at low energies, the parameters have to be hxed by low energy 
quantities like hadron masses, while for high energy perturbative calculations, 
one has to take high energy quantities like jet cross sections. It is not a priori 
clear if one set of hxed parameters describes the ehects of QCD both at low and 
at high energy scales. 

To study this question is one aim of the ALPHA collaboration. To this end, 
one has to dehne the renormalised parameters at low energies on the lattice and 
then evolve them to high energies, where they can be compared to the quan¬ 
tities dehned perturbatively. The technique developed for this purpose is the 
Schrodinger functional scheme , which mainly amounts to dehning the renor- 
malised coupling and masses in a four dimensional box of box size L with special 
boundary conditions, imposing a constant colour electric background held. One 
thus gets a running coupling and running quark masses depending on the scale 
1/L. Having computed these quantities at some low energy scale, they may be 
evolved to higher energies step by step, using the so called step scaling function. 
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In the Schrodinger functional scheme, the renormalised strong coupling is 
constructed by choosing the boundary conditions to be dependent on a parameter 
rj. The renormalised coupling is then dehned as the derivative of an effective 
action with respect to rj. For this reason, a non-vanishing background held is 
required for the renormalisation of the coupling. 

The renormalised quark mass is dehned using the PCAC relation, which con¬ 
nects the axial current and the axial density P“(x) via 

= 2mP°-{x). (1.1) 


A renormalised quark mass can be dehned using appropriate correlation func¬ 
tions containing the renormalised axial current and density. In contrast to the 
continuum case, where the axial current does not need to be renormalised, it does 
get a hnite renormalisation in lattice QCD. This is a special consequence of the 
lattice regularisation, which explicitly breaks the chiral symmetry of the massless 
continuum theory. The renormalisation needed on the lattice is, however, hnite 
and scale independent. The only divergent and scale dependent renormalisation 
constant needed for mass renormalisation is the one of the axial density P°‘{x). 
Thus, the running of the mass is completely described by this renormalisation 
constant. In contrast to the coupling, the mass renormalisation does not refer 
to the boundary gauge helds and may thus be done both with a vanishing or a 
non-vanishing background held. 

In this scheme, the renormalised coupling and masses can be computed non- 
perturbatively by Monte Carlo simulations. Because of the high computational 
costs of simulations in full QCD, hrst results have been obtained in the quenched 
approximation |^, where the fermion determinant is constant, which means 
that the havour number is set to N{ = 0. There have, however, been recent 
results in full QCD with two havours m- A less expensive toy model are so 
called bermions, which is a theory with Af = —2 pSj ^ |30i |3T| . Although this 
is of course no realistic model of the physical world, bermions may be used to 
study structural properties of the theory. The original idea of extrapolating from 
Af = —2 and Af = 0 to Af = 2 did, however, not turn out to be practicable. 

In the high energy regime, one wants to compare the renormalised parameters 
dehned in the Schrodinger functional scheme to those commonly used in this 
energy region, i.e. to renormalised parameters in some perturbative scheme. For 
this purpose, it is necessary to expand the renormalised coupling and masses in 
powers of the bare coupling. For the coupling, this has been done both at 1- 
and 2-loop order 


For the mass, however, there has only been 
a 1-loop calculation with a vanishing background held [^. One of the aims of 
this thesis is to do this calculation with a non-vanishing background held. The 
computation will give the hnite part of the renormalisation constant Zp at 1-loop 
order as well as the 2-loop anomalous dimension. 

Another useful application of perturbation theory on the lattice is the estima¬ 
tion of discretisation errors. For the coupling, the discretisation errors of the step 
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scaling function have been calculated in [Q up to 2-loop order. The 2-loop coef¬ 
ficient does, however, contain the critical quark mass, at which the renormalised 
mass vanishes, at 1-loop order. In [^, the known continuum limit of the critical 


mass was used. In order to get a precise result for the discretisation error, it is, 
however, necessary to use the critical mass at the hnite lattice spacing at which 
the step scaling function is computed. The calculation of the critical quark mass 
and the resulting discretisation errors of the step scaling function of the coupling 
will be done in this theses. 

For the quark mass, several lattice artefacts may be considered. One way to 
estimate the size of the lattice effects is to construct several different unrenor- 
malised masses using different correlation functions containing the axial current 
and density. The difference of these masses should then be a lattice artefact of 
order a (or in the improved theory). For the renormalised quark mass, the 
deviation of the step scaling function from its continuum limit is of interest. This 
deviation has already been calculated with a vanishing background field. If one 
wishes to compute the quark masses in the same runs as the coupling, it is, how¬ 
ever, desirable to have an estimate for the discretisation errors of the step scaling 
function of the running quark mass with a non-vanishing background held. These 
discretisation errors will be computed in this thesis at 1-loop order. 

The thesis is structured as follows: Chapter is an introduction to the basic 
ideas of lattice QCD. Also Symanzik’s 0(a)-improvement is presented here. In 
chapter the Schrodinger functional is dehned, and the renormalisation of the 
coupling constant in this scheme is outlined. The Schrodinger functional is then 
treated perturbatively in chapter including the details of the gauge hxing 
and the expansion of the coupling. Chapter ^ describes how to calculate the 
PCAC mass, serving as a mass later renormalised multiplicatively, on the lattice. 
The critical bare mass, at which the lattice PCAC mass vanishes, is calculated 
at 1-loop order by expanding correlation functions containing the axial current 
and density. In chapter |^, the renormalisation of the mass is studied. The 
renormalisation constant of the axial density and its step scaling function are 
calculated at 1-loop order as well as the discretisation errors. The results of 
all calculations can be found in chapter ^ The thesis closes with a summary 
in chapter ^ Finally, the appendices contain notational conventions, some 
computational techniques used in the calculations, and some numerical results. 







Chapter 2 

QCD on the lattice 


2.1 Lattice gauge theory 

In this chapter, the basic concepts of lattice gauge theory will be introduced. 
Of course, only a short overview can be given here, with an emphasis on ideas 
needed later on in this thesis. More details on quantum helds on the lattice can 
be found in the literature 

Lattice gauge theory is set up on a 4-dimensional hyper-cubic lattice with 
lattice spacing a. In order to bring quantum held theories onto the lattice, one 
has to discretise the action of the corresponding Euclidean continuum theory. 

In the continuum, the gauge helds are represented by the vector potential 
lying in the Lie algebra of the gauge group.Q The calculations in this 
thesis will be done for the case of QCD, where one has the gauge group SU{3), 
but the construction of the theory is applicable for a general gauge group SU (iV). 
On the lattice, the gauge helds are expressed as parallel transporters between the 
lattice points. Let p be a unit vector in direction /i. Then the gauge held on the 
link between the lattice site x and the lattice site x + afi is represented by a link 
variable U{x,^), which is related to the continuum gauge held A^{x) by 

U{x,fj.) (2.1) 

Local gauge transformations are represented by gauge functions ^(a:) living on 
the lattice sites x and acting on the link variables according to 

U (x, jj) —>• U^{x, jj) = Vl{x)U{x^ jj)VL{x -f- (2.2) 

Clearly, U (x, /i) is an element of the gauge group, allowing us to take products 
of link variables on a curve consisting of several links. A special curve giving a 
gauge invariant combination of link variables is a closed loop consisting of four 

^Here and in the following, Greek letters denote Lorentz indices from 0 to 3, while Latin 
indices run from 1 to 3. 


0,0, m, n, n 


5 
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links in the /i-z/ plane, which is called a plaquette p. The product of the link 
variables around this plaquette is then denoted by 

U{p) = U (x, jj.)U{x + a/l, v)U{x + az>, p)~^U (x, (2.3) 


According to Wilson [^, an appropriate discretisation of the gauge field action 
is given by 

= ( 2 . 4 ) 

90 r, 


where the sum is to be taken over all oriented plaquettes. The Wilson action can 
be shown to coincide with the continuum action at leading order in the small a 
expansion. 


2.2 Fermions on the lattice 

2.2.1 The naive fermion action 

In contrast to the gauge fields, the lattice fermions are not situated on the links 
but on the lattice sites themselves. They are Grassmann valued helds and carry 
Dirac, colour, and flavour indices. Thus they differ from the continuum case 
merely by the fact that they are only dehned on discrete lattice sites. 

To set up the theory, one has to define an action for the fermions on the 
lattice. For this purpose, the covariant derivative has to be carried over to the 
lattice. This may be done using the forward derivative, 

Vf.'ipix) = -[U (x, p.)'ip{x + afi) - '^((e)] (2.5) 

Qj 

or the backward derivative, 

V*V'(a:() = - -U{x- ajl, /x)" V(a^ - a/i)] . (2.6) 

CL 

Using the average of these two derivatives as a discretised version of the covariant 
derivative, one gets the naive fermion action 

5'F,naiVe[U,Vz,'0] = '^'ljj{x){D + mo)'ijj{x), (2.7) 

X 

where mo is the bare quark mass. Strictly speaking, mo is a diagonal matrix 
containing the bare masses of the different quark flavours. In this thesis, however, 
only degenerate quark masses will be considered. D is the Dirac operator, 

D = li, (K + K). 


(2.8) 
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where 7 ^ are the Dirac matrices, which can be found in appendix Of course, 
one may add terms that vanish in the continuum limit. While ( p. 8 |) seems to be 
the simplest formulation, it turns out to cause some severe problems outlined in 
the next subsection. For later use, a more complicated expression for the fermion 
action will be needed due to these difficulties. 


2.2.2 Fermion doubling and chiral symmetry 


The main problem of the naive fermion action ( p. 8 |) is the phenomenon of fermion 
doubling. This property is easily seen when writing down the propagator one gets 
from the naive action in momentum space. 


S ‘(p) = mo + - y] sm(p 


(2.9) 


This expression has not only one but 16 zeros in the Brillouin cell in the chiral 
limit mo —> 0 , and these spurious poles of the propagator cannot be ignored, 
since they survive in the continuum limit. To make things even worse, they 
come in pairs with opposite axial charge, thus spoiling the axial anomaly. As a 
consequence of these problems, the naive action is obviously not acceptable as a 
discretised fermion action. 

The occurrence of fermion doublers can also be understood as a consequence 
of the symmetries of the naive action. In the continuum, the flavour symmetry of 
the massless action on the classical level is f/(l)L x D(1 )r x SU{Nf)i^ x S'i7(W)R, 
where Nf is the number of flavours. As a consequence, one has the symmetries 
f/(l)v and S'f/(Af)v, which are valid even in the case of non-vanishing masses, 
and the axial symmetries U (1)a and SU (Af)A- While SU{N{)y and U{l)y remain 
exact symmetries of the quantised theory giving rise to isospin and baryon number 
conservation, the SU{N{)p^ is spontaneously broken, leading to (A^^ —1) Goldstone 
bosons, which explains the relative lightness of the pious. The U{1)a is broken 
in the quantised theory by the axial anomaly. The crucial point here is that the 
naive fermion action (| 2 . 8 |) , instead of x 17(1)^, has a much larger symmetry 

group U (4)l X 17(4)r, where the spurious symmetry transformations exchange the 
corners of the Brillouin cell. 

Several ways have been tried to circumvent the problem of fermion doubling. 
One possibility is to dehne only one spin component per lattice site, thus reduc¬ 
ing the number of fermions to four, which may then be interpreted as different 
flavours. These staggered fermions introduced by Kogut and Susskind M, ^ 


have a major drawback in the fact that gauge interactions break the flavour sym¬ 
metry at hnite lattice spacing. As a consequence, the 16 degrees of freedom on 
the lattice become a mixture of spin and flavour, making the interpretation of 
operators in terms of spin and flavour non trivial. 

Another possibility is to add a dimension hve operator to the fermion action. 
This term will vanish in the continuum limit and can be chosen such that the 
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15 spurious flavours become infinitely heavy. A convenient choice was proposed 
by Wilson |^. Instead of taking the naive Dirac operator ( |2.8|) , the action of 
Wilson fermions is constructed using the Dirac-Wilson operator 


0=5 


v„)-av;v„}. 


( 2 . 10 ) 


Using this action, the spurious fermion states get a mass proportional to 1/a and 
thus decouple in the continuum limit. 

While the problem of fermion doubling is solved by the Wilson action, it in¬ 
troduces a new difficulty. The chiral symmetry of the massless theory is explicitly 


broken by the additional term in (|2.10|) . This is easily seen by stating that 


75 O + 075 ^ 0 . 


( 2 . 11 ) 


This is, of course, a major drawback, because it disables us from treating the¬ 
ories in which chiral eigenstates play a crucial role on the lattice. For these 
purposes, one would like to have a theory free of doublers and preserving chiral 
symmetry. The question whether such a theory does exist is answered by the 
Nielsen-Ninomiya theorem K 


If one uses an action which 


• is translation invariant, 

• has continuous lattice momenta in the range [0, 27r] for L —> 00 , 

• has only local interactions, 

• gives the correct continuum propagator in the continuum limit, 

• and preserves chiral symmetry at finite lattice spacing a. 


then the theory will have doublers. This means that one has the choice either to 
get rid of the doublers or to preserve chiral symmetry, but one can not have both 
at the same time. 

A recent development in this held is the resurrection of the Ginsparg- 
Wilson relation |4^ for a chirally invariant formulation. If the Dirac operator 
satishes the condition 

75D -|- D75 = aD ' j ^ D , (2-12) 

then it has a symmetry that becomes the chiral symmetry in the continuum limit. 
Liischer has presented a detailed analysis how to avoid the implications of the 
Nielsen-Ninomiya theorem using Ginsparg-Wilson fermions 


47 


For our calculations, which do not depend on the chiral symmetry of the the¬ 
ory, it is however simpler to use the Dirac-Wilson operator. The only implications 
of broken chiral symmetry will be a hnite scale independent renormalisation of 
the axial current and an additional renormalisation of the quark mass. 

























2.3 Symanzik’s improvement programme 


9 


2.3 Symanzik’s improvement programme 


The main purpose of lattice gauge theory is to calculate physical quantities on 
the lattice and then extrapolate to the continuum limit. However, the closer one 
gets to the continuum limit, the more expensive are Monte Carlo simulations due 
to critical slowing down. Therefore, it is desirable to make the theory converge 
faster. This is possible due to the fact that the lattice discretisation of the action 
is not unique. One may add irrelevant terms that vanish in the continuum limit. 
This property may be used to cancel 0(a) effects. At small lattice spacing a, one 
gets an effective action 

S-eff = So + aSi + 0(a2) (2.13) 


by Taylor expansion. The principle idea of Symanzik P3|, |2^ was to subtract the 
order a term, thus making the theory converge at a rate proportional to a^. 

The pure gauge action can be shown to reach its continuum limit at O(a^) 
without any additional terms, but the quark action does need improvement. For 
on-shell quantities, the 0(a) contribution can be cancelled by the Sheikholeslami- 
Wohlert term 1081 


V', V'] = '0(a;)^T)v'0(a;), 


(2.14) 


with 




The Dirac structure of the Sheikholeslami-Wohlert term is given by 


(2.15) 


(2.16) 


and is the lattice field tensor given by 

F^,y{x) = ^ {Q^u{x) - Qu^^{x)} , (2.17) 

with 

Qfiv{x) = |l7(x,/i)t/(x + a/t, i/)t/(x + az>, ya)“^17(x, 1 /)“^ 

+U{x, u)U{x + ah' — afi, {x — afi, v)~^U{x — afi, jj) 

+U{x — afi, {x — aO — afi, v)~^U{x — au — afi, fi)U{x — aO, u) 

+U{x — ah', i>)~^U{x — at), fi)U{x — aP + afi, z/)17(x,/i)“^|. (2.18) 

The Sheikholeslami-Wohlert term may be visualised by the characteristic shape 
in figure |2.1| . Therefore, it is often called the clover term. The improvement 
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V 




Figure 2.1: The Sheikholeslami-Wohlert term 


coefficient Csw may be computed in perturbation theory. Up to 1-loop order, the 
result originally obtained by Wohlert m is 


C.W = 1 + 0.26590(7)g„' + 0(gj), 


(2.19) 


independent of Ni. Values obtained later using the Schrodinger functional ^ 
are completely compatible with this result. For iVf = 0 and N{ = 2, Csw has also 
been computed non-perturbatively |^, For 0 < (70 < 1, the results can be 
represented in good approximation by 


Csw(5'o) lArf =0 

Csw {do) I Ni=2 


1 - 0.656^g - 0.152c/^ - 0.05%^ 

1 - 0.922C/2 ’ 

1 - 0.45%g - 0.175c/^ + 0.012C/6 + 0.045c/^ 
1 - 0.720172 


( 2 . 20 ) 

( 2 . 21 ) 


For the case of Nf = —2 (“bermions”), one gets good results by linear extrapola¬ 
tion of the results for N{ = 2 and iVf = 0 [^. 

However, the improvement of the action is not sufficient to improve the ex¬ 
pectation value of every composite operator. The operator has to be improved 
in addition. To achieve this improvement, one has to expand the operator in 
powers of a and then remove the 0(a) contribution. To this end, one has to 
write down a basis of operators with the correct dimensions and the symmetries 
required and subtract them with appropriate improvement coefficients. Often, 
the number of operators can be reduced by using the equations of motion. An 
important example is the improvement of the axial current |5^, which will play 
an important role in later chapters. 
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2.4 Renormalised parameters 

It is well known that, in order to get finite results, quantum field theories have 
to be renormalised, leading to a redefinition of the parameters of the theory, like 
couplings and fermion masses. The renormalisation is done by first introducing 
a regularisation, then doing the renormalisation cancelling divergences, and fi¬ 
nally removing the cutoff again. This situation does also occur in lattice QCD, 
where the cutoff is given by the inverse lattice spacing a~^. In the following, the 
renormalisation of the coupling and the masses will be briefly summarised. 


2.4.1 The renormalised coupling 

Since, in the continuum, the coupling can only be calculated perturbatively, one 
has to apply a perturbative renormalisation scheme to compute the renormalised 
coupling g. Such schemes are the MOM-scheme [^, the MS-scheme [^, and 


the MS-scheme of dimensional regularisation. The last two schemes differ 
by the subtraction of a constant in the MS-scheme in addition to the divergence, 
while in the MS-scheme, only the divergence itself is removed. This variety of 
renormalisation schemes is in contrast to the case of QED, where one naturally 
chooses the scheme such that the charges and masses of the leptons get their 
physical values. This is however not possible in QCD, because free quarks cannot 
be observed. 

On the lattice, one can define a renormalised coupling non-perturbatively 


One possibility is to consider a pair of a static quark and a static antiquark 
separated by the distance r. The force F{r) acting between these quarks is the 
derivative of the static potential V (r), which can be computed from Wilson loops. 
Defining a = g'^/ATi, a physical coupling is given by 


Of 


( 2 . 22 ) 


where q = and C-p is the eigenvalue of the quadratic Casimir operator of the 
gauge group in the fundamental representation. For SU{N), it is 


Cf = 




2N 


(2.23) 


For QCD, one thus has Cp = 4/3. 

What all couplings computed in different renormalisation schemes have in 


common is the fact that they depend on a renormalisation scale, like q in (|2.22|) . 
The running of the coupling is described by the Callan-Symanzik equation 
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where the /3-function has an asymptotic expansion 

Pid) ~9^ [^0 + + 0 (^^)] . ( 2 . 25 ) 

In general, the expansion coefficients will be scheme dependent. The first two, 
however, are the same in any two mass-independent renormalisation schemes. In 
these schemes, the renormalisation conditions are imposed at zero quark mass, 
thus avoiding an implicit dependence of the renormalised coupling and fields on 
the quark mass. Then ho and hi are given by 

^ ( 4 ^ “ 3 ^') ’ 

For energies high enough to make perturbation theory possible and for sufficiently 
small flavour number iVf, one has asymptotic freedom, which means that the 
coupling approaches zero in the limit of infinite energy. The asymptotic solution 
is 

5^00 1 _ fei ln[ln(gVA^)] / {ln[ln(gVA^)]}^ \ 

^ ^ &oln(g7A2) hlMqyA^)]^ \ [ln(g7A2)]3 ^ 

where A is a scheme dependent integration constant. In the high energy regime, 
it may be used to relate different renormalisation schemes to each other. 


2.4.2 The renormalised quark masses 

Like the coupling, also the quark masses have to be renormalised. In continuum 
perturbation theory, their renormalisation is incorporated in the MOM-, MS-, or 
MS-scheme. On the lattice, a good choice is the hadronic scheme, in which the 
bare quark masses are eliminated in favour of physical hadron masses. First, one 
chooses certain values for the bare coupling go and the bare masses am^, where the 
index f labels the different quark flavours u,d,s,c,b. Neglecting isospin breaking, 
one may assume the light quarks to be degenerate and define mj, = mj) = mj). 
Next, one calculates the masses of five different hadrons H containing quarks of 
all flavours that are to be renormalised, for example H=p,7r,K,D,B. The hadron 
masses will, of course, depend on the bare parameters, 

amn = amii{go, amj,, ami, arn^). (2.29) 

To renormalise the theory, one first sets the proton mass rUp to its experimental 
value determining the lattice spacing by 

amp 


a = 


(2.30) 
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Next, one has to choose the parameters amg such that the hadrons get the masses 
known from experiment. Equivalently, one can hx the bare quark masses at a 
given value of go from the condition 


amn am 


exp 

H 


amr 


am 


exp 


(2,31) 


with H = 7r,K,D,B. The bare coupling then determines the lattice spacing via 2.30 


Once having renormalised the masses this way and the coupling according to p.22 
the theory is completely dehned in terms of physical observables. 

Like the renormalised coupling, the renormalised masses are scale dependent. 
To describe their running, one may use the r-function. 


dm _ 

= ngW 


(2.32) 


It has an asymptotic expansion 


r{g) -g^ [do + g^di + 0{g^)] , 


(2.33) 


where all expansion coefficients are scheme dependent except the 1-loop anoma¬ 
lous dimension do- It is given by 


do — 


(dvr) 


(2.34) 


The 2-loop anomalous dimension di is known in the MS scheme. For the gauge 
group SU(N) it is given by 


^ (47r)4 \ 6 2 3 ' 


(2.35) 


As an analogy to the A-parameter, one may introduce a renormalisation group 
invariant quark mass M, dehned by 


M = lim m{2bog‘^) 


q—>^oc 


(2.36) 


It can be shown to be independent of the renormalisation scheme. It can be used 
to obtain the running mass in any scheme by hrst computing it in one scheme 
and then inserting the proper jS- and r-functions of the other schemes into the 
renormalisation group equations |58[ . 

A and M are fundamental parameters of QCD, which are not hxed by the 
theory, but have to be determined from experimental results. 
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2.4.3 Finite renormalisations 


Having renormalised the coupling and the masses, the running parameters at the 
momentum scale q can be obtained by integrating (p.24|) and ( 2.32 ) with the 
boundary conditions 




(2.37) 


Any two mass-independent renormalisation schemes are related by a hnite 
renormalisation of the parameters, 


/i' 

o' 

A 

(2.38) 

9r 

- 9R\/xgi9R), 

(2.39) 


= m^Xm(^R)- 

(2.40) 


The hnite renormalisation constants Xg and Xm may be expanded in perturbation 
theory, 

OO 

= + ( 2 . 41 ) 

k=l 

With this expansion, one hnds that the 2-loop anomalous dimensions in both 
schemes are related by 


d[ = di + 2box^r^^ - doXg^■ (2-42) 

The crucial point of this equation is that and are 1-loop coefficients. This 

means that, knowing the 2-loop anomalous dimension in one scheme, one may 
compute it in a different scheme by a 1-loop calculation. This will be used later 
in order to get the 2-loop anomalous dimension in the Schrodinger functional 
scheme from the 2-loop anomalous dimension in the MS-scheme. 






Chapter 3 

The Schrodinger functional 


3.1 Motivation 


In the last chapter, it was seen that one has two kinds of renormalised couplings: 
The perturbatively dehned couplings like extracted from high energy exper¬ 
iments, and Oqg, dehned non-perturbatively at the scale of hadron physics. Now, 
an obvious task for lattice QCD is to connect these couplings. This amounts to 
calculating the coupling non-perturbatively at low energy scales and then evolv¬ 
ing it to high energies, where one may compare it with high energy experiments 
like the determination of jet cross sections. Thus one should be able to verify 
that the hadron spectrum and the properties of jets really are described by the 
same theory. 

To do this calculation, one has to fulhl several requirements at the same time. 

• One has to calculate agg{q) at energy scales g ~ 10 GeV or higher, in order 
to make the connection to other schemes with sufficiently small perturbative 
errors. 

• One has to keep the renormalisation scale q sufficiently far from the cut¬ 
off a~^ in order to keep the discretisation errors small. Otherwise, a safe 
extrapolation to the continuum limit might not be possible. 

• One has to keep the box size L large compared to the conhnement scale in 
order to avoid large hnite size effects. 


These criteria can be summarised by 


r 1 1 

L > TT-TT^ > - 

0.4 GeV q 


10 GeV 


> a. 


( 3 , 1 ) 


This means that one has to perform a Monte Garlo simulation on a four dimen¬ 
sional lattice with a size of L/a S> 25. Even with modern computer technique, 
lattices of this size are far from being accessible. 
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i^max = O(0.5fm) : HS SF(g = l/L^^x) 

i 

SF(g = 2/L^ax) 
i 


i 

SF(g = 


PT: i 

PT 

jet - physics <— Aqcd, M 

Figure 3.1; Strategy for evolving the renormalised coupling to large energy scales 


However, there is a solution to this problem, 
by identifying the two physical scales ||2^, 


The difficulties can be avoided 



(3.2) 


This means one takes a hnite size effect as the physical observable. The strategy 
one may use now is presented in figure |3.1| . First, the hadronic scheme HS 
can be related to a finite volume scheme denoted by SF at a low energy scale 
q = 1/Lmax- As a finite volume scheme, we will use the Schrodinger functional, 
which will be defined in the next section. Within this scheme, one may now 
compute the evolution up to the desired scale q = 2"^/Ljnax step by step. With 
this recursive procedure, one can thus avoid using quantities defined at energy 
scales very wide apart. At an energy scale large enough for perturbation theory, 
one can then compute the A-parameter and the renormalisation group invariant 
quark masses in order to make the connection to renormalisation schemes defined 
perturbatively. 


3.2 Definition of the Schrodinger functional 

The purpose of this section is to define the Schrodinger functional and introduce 
the basic notation. This will be done in the continuum. In the next section, these 
definitions will be carried over to the lattice. 
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3.2.1 Formal definition 

We start with the definition of the Schrodinger functional in pure gauge theory 
and later add the fermions to this framework. 

The Schrddinger functional is basically defined by the Hamiltonian evolution 
of the gauge fields. To this end, one has to specify the theory at xq = 0 (or any 
other fixed time) and write down the Hamiltonian of the theory. 

Since the main purpose of the Schrodinger functional method is to study the 
scaling properties of QCD at a finite volume, the theory is set up in an L x L x L 
box where L serves as the scale, and one assumes periodic boundary conditions for 
the gauge fields. The gauge fields are represented by vector potentials with values 
in the Lie algebra of SU{N). The calculations in this thesis will be done for the 
case of SU{3), but the definition is applicable for arbitrary N. In the temporal 
gauge, one is thus left with Lie algebra valued vector components ^^(x) which 
have to be periodic.^] Under gauge transformations A(x), they transform like 

Hfc(x) ^ Hfc(x) = A(x)Hfc(x)A(x)“^ + A(x)4A(x)"b (3.3) 

In order to preserve the periodicity of the gauge fields, only periodic gauge func¬ 
tions A(x) may be allowed, which can be regarded as continuous Lie group valued 
functions on a 3-dimensional torus. These functions are topologically non-trivial 
and fall into disconnected classes that are distinguished by an integer winding 
number. If one wrote down the Schrodinger functional in the functional inte¬ 
gral representation, where one has to integrate over all gauge configurations, one 
would have to sum over all topological classes. However, since this difficulty does 
not occur on the lattice, it will not be studied in detail here. 

Now, the quantum mechanical states may be defined as wave functionals '(/'[A] 
acting on the gauge fields. Then the quantity 

im = j D\A]^\A]-x\A] (3.4) 

defines a scalar product on these states. The measure D[A\ is meant to be 

D\A] = n (3.5) 

where A^(x) are the components of Afc(x) in a basis of the Lie algebra su{n). 
The basis used in this thesis is explained in appendix E However, not all states 
defined this way are physical. Physical states have to be gauge invariant, which 
means they have to satisfy 

= V'iA]. (3.6) 


^Here and in the following, bold letters like x denote vectors in the three dimensional space. 
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Obviously, these physical states form a subspace of the space of states. Any state 
can be projected onto this subspace by a projector P, given by 

¥^/J[A] = j D[A]^/j[A% (3.7) 

where the measure 71[A] is dehned by 

D[A] = JJdA(x). (3.8) 

X 


The canonically conjugate held of the gauge held is the colour electric held 


(3.9) 


It is part of the colour held tensor, with the magnetic components given by 

F“(x) = d,A?(x) - d,A“(x) + r^At(x)A?(x), (3.10) 

where are the structure constants of SU(N). Using this tensor, one can write 
down the Hamilton operator H, 




(3,11) 


Now, for any smooth gauge held C^, a state \C) may be dehned such that 


{C\i,) = m (3,12) 

holds for all wave functionals This state is not necessarily physical, but of 

course it can be made so by projecting it onto the physical subspace using the 
projector P. The Schrodinger functional Z[C'^C] is then dehned by 

Z[C\C] = {C'\e-^^^C). (3.13) 


The projector P makes sure that the Schrodinger functional is invariant under 
gauge transformations of C and C", since only gauge invariant intermediate states 
contribute. 

An alternative way to write down the Schrodinger functional is to set up the 
theory on a 4-dimensional box L x L x L x T and express the time evolution 
operator by a functional integral, using the action 

/ d^a:tr{F^^F^^}. (3.14) 

^9o Jo 

The helds C and C then dehne the boundary conditions on the gauge held at 
xo = 0 and xq = T. This situation may be depicted by hgure On the lattice, 
this representation will be chosen. 
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LxLxL 


Figure 3.2: The space time box of the Schrodinger functional 


3.2.2 Fermions in the Schrodinger functional 

Though mostly used on a lattice, also the fermion part of the Schrodinger func¬ 
tional may formally be defined in the continuum [^, |^. The easiest way to do 
this is in the functional integral representation with the action 




dxo / d^xi/j{x) {'yf.Df, mo) -ipix) 
Jo 

r d^x [il}{x)P-7p{x)]^^^^ 
f d^x [4J{x)P+^|J{x)] , 


(3.15) 


where = d^ + is the covariant derivative of the fermion fields. The projec¬ 
tors P+ and P_ are defined by 


= 2 (1 ± 7o) ■ 


(3.16) 


In time direction, Dirichlet boundary conditions are imposed for the quark fields 


P+V'(a:) 1^0=0 = p(x), 
and for the antiquark fields 

^{x)P-\xo =0 = p(x). 


P-'lpix)\^Q=T = p'(x). 

(3.17) 

V'(a:)P+|xo=T = p'(x). 

(3.18) 
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In space direction, the qnark fields are chosen to be periodic up to a phase factor, 

+ Lk) = '4){x + Lk) = e~'‘^'4){x). (3.19) 

For the time being, 9 will be left arbitrary. Later, several values for 9 will be 
considered. 


3.3 Lattice formulation 


3.3.1 The Schrodinger functional action on the lattice 

Since the Schrodinger functional is characterised by its boundary conditions, these 
boundary conditions have to be carried over to the lattice. This means that the 
link variables have to take fixed boundary values in the Xq = 0 and Xq = T planes, 

U{x, k), U{x, k)\,,=T = k). (3.20) 


In order to make contact to the continuum definition of the Schrodinger func¬ 
tional, one has to relate the lattice boundary fields W and W to the continuum 
boundary fields C and C. This relationship is naturally established by recalling 
that the link variables are parallel transporters for colour vector fields. So it 
seems natural to define 


W (x, k) =V exp 


dt Ck{'^ + ak — tak) 


(3,21) 


and analogously for W using the continuum field C instead of C. Here, V 
denotes the path ordering along the link. 

The gauge field action in the Schrodinger functional is given by the Wilson 
action 

P-22) 

9o p 

only modified by the weight factors w{p). In the unimproved theory, one has 
w{p) = 1/2 for the spatial plaquettes at Xq = 0 and Xq = T and w{p) = 1 in 
all other cases. In the improved theory, these weight factors will be modified at 
the boundaries. Here again, the sum in (|3.22|) is to be taken over all oriented 
plaquettes. 

For the quarks, we take the usual Dirac-Wilson action. 




with the Dirac-Wilson operator (|2.1CI|) 


^ + V^) - aV/V4 . 


(3.23) 


(3.24) 
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The Schrodinger functional on the lattice is then given by 

= j (3.25) 


where the action S is the sum of the gauge field action and the fermionic action, 
S = Sq + Sp. The expectation value of any product O of fields is now given by 


{O) = [ D[U]D[2jj]D[ijj]Oe-^^^’^’^A . (3.26) 

V J p'=p'=p=p=0 

Apart from the quark and gauge fields, the operator O may also contain the 
boundary fields 


C(x) 

ax) = 


_ 6 _ 

5p(x) 

S 

(5p(x) 


C'(x) 

a(x) 


5 


(5p'(x)’ 

6 


(5p'(x) ■ 


(3.27) 

(3.28) 


These derivatives have to be understood in the sense that they act on the Boltz¬ 
mann factor in (|3.26|) . They will be used to define correlation functions which 
play an important role in later chapters. 


3.3.2 The background field 


Unless they are zero, the boundary fields C and C" obviously impose a background 
field which is given by a solution of the field equations with the correct boundary 
values. It will be explained later that a non-zero background field is inevitably 
needed to define the renormalised coupling in the Schrodinger functional. Since 
only small lattices are accessible for numerical simulations, one cannot make the 
lattice spacing arbitrarily small. Thus one has to deal with considerable cutoff 
effects, which requires a strategy to keep the discretisation errors as small as 
possible. In particular, one has to look for a background field leading to minimal 
cutoff effects. 

It can be shown that spatially constant Abelian boundary fields give good 
results [Effl. So one may choose 


^ f (j)ki 0 0 \ 

Ck = j\ 0 0 , (3.29) 

VO 0 / 


and analogously for C. In the following, and C'f, will be chosen to be inde¬ 
pendent of k, i.e. Cl = C 2 = C 3 and C[ = C 2 = Cg. An obvious solution of the 
field equations with these boundary conditions is 


Bo = 0, Bk = [xoC^ + {T - xo)Ck] /T. 


(3.30) 
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Gkl — 0, 

(3.31) 

- Ct]/T. 

(3.32) 


This solution is a constant colour electric field, which can easily be seen from the 
corresponding field tensor 


with 


On the lattice, the background field must of course be given in terms of the link 
variables. A good candidate for the lattice background field is 

l/(x,/i) = exp{ai?^(a;)}, (3.33) 

which, after inserting ( 3.30|) , gives a field that only in the spatial directions is 
different from 1, 

y(x,0) = l, Vix,k) = Vixo), (3.34) 

with 

V{xq) = exp{i[£xQ — iC]} . (3.35) 

This field trivially fulfils the boundary conditions and it can easily be shown to 
be a solution of the field equations. Yet, this field is still only a candidate for 
the background field. It has to be shown that V really is a configuration of least 
action and that the same action is only given by fields that are gauge equivalent 
to V. According to |^, this is the case for the gauge group SU{N) if the (j)ka 
and the lie in the so called fundamental domain 

N 

4>kl < 4>k2 < . . . < \4>ka “ 0fc/3| < 27r, ''^^(pka = 0, (3.36) 

a=l 

and the relation 

TL/a^ >{N- l)7rVax{l, N/16} (3.37) 

holds. This relation is only needed in the proof of the theorem for technical 
reasons, it is of no physical significance. For N < 3 and the lattices of interest, 
it is unimportant. 

A convenient choice of boundary fields obeying the conditions mentioned 
above is 

0 




0 

0 


r7(-i + i/) 


0 

0 


0 


-r/(i + z/) + f 


(3.38) 


for the boundary field at Xq = 0 and 

-T] 


Ck = T 
^ L 


TT 0 0 

0 v{\ + ^) + l 0 

0 0 + 


(3.39) 


■ 2 ' 2 , 

for the field at Xq = T. Here, we will choose i/ = 0. Later, also rj will be set 
to zero, but first it is needed as a variable in order to define the renormalised 
coupling in section [O. 
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3.4 0(a) improvement of the Schrodinger func¬ 

tional 

In the Schrodinger functional, one has two kinds of counter-terms for 0(a) im¬ 
provement. The hrst one is the volume term, improving the helds in the interior 
of the lattice. The second kind are the boundary counter-terms arising from the 
Schrodinger functional boundary conditions in the time direction. The improved 
action thus becomes a sum 

SimAu, A = S[U, ij, A + 5SAU, A + sScAA + ss^Au, A, (3.4o) 

where the volume term SSy is the Sheikholeslami-Wohlert term familiar from 


section 2.3 


On the boundaries at Xq = 0 and Xq = T, composite helds of dimension 4 
have to be subtracted. There are contributions of order a both from the gauge 
helds and the quark helds 1^^. The gauge held contribution is cancelled by the 
counter-term 

■S'Sg.i.IC'I = 5^(c.-l)^tr{l-f/(p,)} 


Ps 


--2(ci-l)5^tr{l 


U(Pt)}, 


(3.41) 


where the sum runs over all oriented plaquettes at the boundaries, both time- 
like {pt) and space-like (ps). This boundary improvement term simply amounts 
to changing the weights w{p) in (|3.22|) to w{pt) = Ct and ta(ps) = Cs/2. The 
boundary improvement coefficients have perturbative expansions 


Ct = l + cJ^^Po + c?^^o + 0(^o)i 

Cs = l + cAgl + cAgo + 0{g^). 


(3.42) 

(3.43) 


The 1-loop coefficient for the time-like plaquettes can be decomposed into a part 
independent of the havour number and a term linear in iVf, 


,0) - „(■.«) + 


where the coefficients are given by |^, ^ 


= -0.08900(5), 
= 0.0191410(1). 


^(1,1) 

Ct 

The 2-loop coefficient has a decomposition 

cf) = cf’n cf 


(3.44) 


(3.45) 

(3.46) 


(3.47) 
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with [^, H 


.(2,0) 

t 

= -0.0294(3), 

(3.48) 

.(2,1) 

't 

= 0.002(1), 

(3.49) 

.(2,2) 

't 

= 0.0000(1). 

(3.50) 


With the special choice of spatially constant Abelian boundary helds, the im¬ 
provement term for the space-like plaquettes vanishes and Cg may be disregarded. 

For the boundary counter-terms depending on the quark helds, all possible 
composite helds of dimension 4 and how most of them are eliminated are presented 
It can be shown that one is left with four terms which, ensuring time 


m 


reversal invariance, are given by 


iSF,b|f/y.V>] = a‘‘5^{(c.-l)|0.(x) + 6;(x)| 



(ct - l)[C>t(x) + C>((x)]|, 

(3.51) 

with the helds 



d(x) 

= ^P(x) 7 fc(Vfc +Vfc)p(x), 

(3.52) 

(^s(x) 

= ^P'(x) 7 fc(Vfc +Vfc)p'(x), 

(3.53) 

(^t(x) 


(3.54) 


= |^(l/)^-Vo'0(2/)+'0(2/)VoP+V'(l/)}^ 

(3.55) 


Setting the boundary quark helds to zero, one is left only with C>t(x) and C>((x). 
The volume term and the quark boundary term may then be rewritten as a 
change in the Dirac-Wilson operator. 


6D = dDy + 5T)b, 


with the volume term 


SDv'ip^x) 


i 




(3.56) 


(3.57) 


and the boundary term 


6D^'tp{x) = 'tp{x)-U{x - ^P+ii{x-a^) 


+ 5; 


XQ^T—a 


'ip{x) — U{x, 0)P-'il}{x -I- aO) 


(3.58) 
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The improvement coefficient Ct can be computed in perturbation theory ^ 
and is, up to 1-loop order, given by 


Ct = l-OM795{2)gl + 0{g^o). 


(3.59) 


Having improved the action, the expectation value of any improved field poly¬ 
nomial O converges to its continuum limit at a rate proportional to a^. First 
integrating over the quark fields, this expectation value may be decomposed into 
a quark field average and a gauge field average. 


(O) = {|01f)g . (3.60) 

where [.. .]f is the quark field average at a given gauge field configuration U and 
(.. .)g denotes the gauge field average with a probability density proportional to 

det{D + 6D + mo) exp {—S'g)?/] — (5*S'G,b} • (3.61) 

The quark field average [0]y can be defined by a generating functional, 

2F[p',p']p,p]V,r];U] = j T>['0]T>['0]exp|-^F,impr[f^,V','0] 

W [ij{x)r]{x) + f]{x)'il;{x)] |, (3.62) 


where S'Fjmpr is the improved fermion action and p{x) and fj^x), 0 < xo < T, are 
source fields for the quark and antiquark fields in the interior of the space-time 
box. Substituting 


'ip{x) 


6 


6p{x )’ 

one gets the quark field average by 


V'(x) 


6p{x )’ 



(3.63) 


(3.64) 


This decomposition into quark and gauge field averages will be used in later 
calculations. 


3.5 The renormalised coupling in the Schrodin¬ 
ger functional 

3.5.1 Definition of the coupling 

Having introduced the Schrodinger functional, we can now turn to defining the 
renormalised coupling in this framework. In order to make practical use of this 
method, certain criteria have to be met by the coupling. 
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Since a central point in the Schrodinger functional method is the wish to 
compute a renormalised coupling non-perturbatively, the coupling should be easy 
to measure in Monte Carlo simulations. Furthermore, the results obtained in 
these simulations should be easy to extrapolate to the continuum limit without 
too large uncertainties. For this reason, discretisation errors should be small. 
And, last but not least, the coupling should have an easy perturbative expansion, 
such that the beta function can be computed to sufficient order and a matching 
can be done to perturbatively dehned renormalisation schemes in the high energy 
regime. 

Considering all these criteria leads to a dehnition based on the effective action 

T[V] =-\nZ[C\C]. (3.65) 


Since the background held V depends on the parameter rj in the boundary gauge 
helds ( p.38|) and (|3.39|) , the effective action will be a function of r], too. So one 
may take the derivative 

r'V] = Arp/]. (3.66) 

and expand it in powers of the bare coupling qq, 


r'lr] = + r;[i/] + gjr' [i/] + 0(g„"). (3.67) 

F' can be shown to be hnite when it is expressed in terms of a renormalised 
coupling like g^. Properly normalised, it may itself serve as a renormalised 
coupling. Thus, the Schrodinger functional coupling is dehned as 


f(L) 




Fp/] 


r )=0 


(3.68) 


The normalisation factor Fg makes sure that the coupling dehned this way coin¬ 
cides with the bare coupling at tree level of perturbation theory. This coupling 
satishes all requirements mentioned above. As the expectation value 


r'M = 



(3.69) 


it is easily accessible to Monte Carlo simulations. The discretisation errors turn 
out to be tolerable. And it is relatively easy to expand in perturbation theory. 
The problem of the perturbative expansion of the Schrodinger functional will be 
addressed in chapter ||. 


3.5.2 The step scaling function and its lattice artefacts 

By dehnition, the Schrodinger functional coupling dehned in this chapter runs 
with the box size L. In order to study its scaling behaviour, one needs a recipe 
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for evolving the coupling from one length scale to another, which means that 
some kind of discrete /3-function is required. An appropriate function in the 
Schrodinger functional is the step scaling function. First of all, one has to start 
from a given coupling u = g‘^{L). Then the length scale is changed by a factor 2. 
(In principle, one may choose a different scale factor, but 2 is commonly used for 
obvious reasons. If one took a larger factor, one would have to perform simula¬ 
tions on larger lattices.) Then, keeping the bare coupling fixed, the coupling at 
the length scale 2L has a value g^{2L) = u' . The step scaling function a is then 
dehned as 

(t{u) = u' . (3.70) 

Obviously, a is an integral over the /3-function. In fact, it has a perturbative 
expansion 

a{u) = u + 2ho ln(2) + ... (3-71) 


However, this is only true in the continuum limit. On a lattice with finite lattice 
spacing a, the step scaling function will depend on the lattice resolution a/L. So, 
instead of (|3.70|) , one has to define 


S(m, a/L) = u'. 


(3.72) 


The continuum step scaling function a is the continuum limit of S, i.e. (t{u) = 
lim.a/L^o'Liu^a/L). In practice, the step scaling function is obtained by simula¬ 
tions on pairs of lattices to determine 'Liu, a/L). One then gets aiu) by extra¬ 
polating these results to the continuum limit. The deviation of the step scaling 
function from its continuum limit. 


Liu,a L) - aiu) 

diu,a L) = --- 

aiu) 


(3.73) 


may be expanded in perturbation theory to give an idea of the size of the cutoff 
effects. Details of the perturbative expansion will be presented in chapter 
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Chapter 4 


Perturbation theory in the 
Schrodinger functional 

4.1 The gauge fixed action 

4.1.1 Preliminaries 

It is a well known fact that, in order to set up perturbation theory, one has to hx 
the gauge. Otherwise, zero modes would occur, since the minimum of the action 
is degenerate in the directions of gauge transformations. 

The Schrodinger functional on the lattice is invariant under all gauge transfor¬ 
mations that leave the boundary helds unchanged. This condition is only satished 
by gauge functions 0(a;) that are constant and diagonal on the boundaries Xq = 0 
and Xq = T. These gauge functions form a group Q. However, not the whole 
group Q needs to be hxed. Those gauge functions f2(x) which are constant and 
diagonal not only on the boundaries but on the whole lattice act trivially on the 
background held. They form a subgroup of Q isomorphic to the Cartan subgroup 
Ctv of SU (A^), which may be factored out and then survives as a global symmetry 
of the theory. The gauge hxing has thus to be done on the group 

Q = Q/Cn. (4.1) 

Q may be identihed with the subgroup of transformations VL & Q that are equal 
to a hxed constant diagonal matrix at xq = T. Here we choose f2(a:)|3;o=r = 1- 
The Lie algebra Cg of Q consists of helds ui with ui^x) G su{N). The gauge 
functions H may be parametrised as 

H(a;) = exp{-5fo‘^(a:)} = 1 - goUj{x) + 0{gl). (4.2) 

Obviously, a; must obey the boundary conditions 

a;(0,x) = fi:, ci;(T, x) = 0, (4.3) 
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where k is constant and diagonal. 

Analogously, one may expand the gauge held U around the background held. 
Let Ti. be the set of gauge helds satisfying the Schrodinger functional boundary 
conditions. In a neighbourhood of the background held, any gauge held U may 
be parametrised by 

U{x,^) = exp{goaqf,{x)}V{x,fi) 

= !^l+goaq^{x) + ^gQa^q^{xf+ 0{gQ)"^V{x,ix), (4.4) 

with q^{x) G su{N). In perturbation theory, the vector helds q become the 
gluons. They form a linear space C'h and obey the boundary conditions 


gfc(0,x) = 0, qk{T,x) = 0. (4.5) 

The time components qo{x) are unconstrained and dehned for 0 < Xq < T. We 
will however get boundary conditions for go later by extending the lattice. 

For later use, it is useful to dehne a scalar product on Cg and For two 
vectors g and r in an obvious choice is 

(g, r) = -2a^ ^ tr{g^(a;)r^(a:)}. (4.6) 

X,ll 

The scalar product on Cg is dehned similarly. 

For the dehnition of a suitable gauge hxing function, we will need the covariant 
forward and backward derivatives on the lattice, 

D^f{x) = -[V{x,fx)f{x +ail)V{x,fx)-^ - f{x)] , (4.7) 

T>*/(x) = ^[f{x)-V{x-afi,n)~^f{x-afi)V{x-ail,n)]. (4.8) 

Using the forward covariant derivative, one may dehne a linear operator d \ Cg ^ 
C-H by 

{du)^{x) = D^u{x). (4.9) 

Another operator that will be used later is d*. It is dehned as minus the adioint 
of d, 

{d*q,uj) = —{q, du) (4-10) 

for all g G C-}^ and u) & Cg. Starting from this dehnition, d* is explicitly given by 

f if 0 < Xq < T, 

(d*g)„g(a;) = <j (aV^^) Ey <?o(0, y)a /3 if a = (3 and Xq = 0, (4.11) 

I 0 otherwise. 
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This rather complicated formulation may be simplified to 




(4.12) 


for all Xq by extending the field q to Xq = —1 and Xq = T and imposing ap¬ 
propriate boundary conditions on the time component go- The new components 
are physically irrelevant and only used for technical reasons. Due to the special 
choice of background field used here, a complication does occur in this procedure. 
If one used, for example, the self dual background field considered in the con¬ 
tinuum formulation in 0 , one would get Neumann boundary conditions for all 
go- This is also the case with the Abelian background field used here, except for 
one case. If go has got a part which is diagonal and spatially constant, this part 
satisfies Dirichlet boundary conditions at Xq = —1. For this reason, the bound¬ 
ary conditions are easier to define in momentum space, where spatially constant 
fields g simply become gluons of zero momentum. So we switch to momentum 
space by taking the Fourier transform 


go(p,xo) = 5^6 *P’‘go(x), 

X 

which may be decomposed in a basis of the Lie algebra of SU{N), 


(4.13) 


<Ir>(p,Xo) = qS(p,Xo)I‘‘ 


(4.14) 


The special choice of the basis /“ for the case of SU{3) used in later calculations 
can be found in appendix 

Now the boundary conditions are given by 


^o?o(PiO) 

d*oro{p,T) 


0 if /“ G Cn and p = 0, 
0 else, 

0 . 


(4.15) 


4.1.2 The gauge fixing procedure 


In the following, we will use the gauge fixing procedure outlined in [^. To this 


end, we will need a gauge fixing function F{U), which is conveniently chosen to 
be a linear mapping from Ti to Cg. This function is defined in a neighbourhood 
M of the background field V and has to fulfil several conditions. First of all, it 
has to vanish on the background field. 


F{V) = 0. (4.16) 

The first order variation of F under a gauge transformation generated by a; is a 
linear operator L{U) ; Cg Cg, 


L{U)u = 5^F{U). 


(4.17) 
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For a gauge fixing function, the determinant of L{U) is required not to vanish. 

A suitable function that fulfils all conditions and is relatively easy to handle 

is 

F{U) = d*q. (4.18) 

Let now f{U) be a function that is non-zero only in the neighbourhood A/”. It 

can be shown that 

I D[U]f{U) = k I D[U]f{U)5{F(U))demu)). (4,19) 

This relation stays valid if one replaces the delta function by 6{F{U) — Z), where 
Z is an element of Cg with (Z, Z) < e, where e is chosen to be appropriately 
small. Since the left hand side of ( |4.19|) is independent of Z, one can integrate 
over Z with a Gaussian measure, resulting in the gauge fixing term 

Ssf|V9l = L(cC,V,) (4.20) 


in the action. 

In the usual way, the determinant in (|4.19|) may be rewritten in terms of Lie 
algebra valued Grassmann variables, the Faddeev-Popov ghost fields c and c with 
the action 

SFp[V,q,c,c] =-{c,d*6cq). (4.21) 


Here, 6cq denotes the first order variation of q under the gauge transformation 
generated by c. Expanding it to order yields 


6cq^, = D^c + goAdqf^c 


^goaAdq^ + ^{goaAdq^f + 




(4.22) 


with no sum over /r implied. 

Also he ghost fields c and c have to satisfy certain boundary conditions. Like q, 
also the ghost field c may be Fourier transformed and decomposed in the basis /“ 
of the Lie algebra of SU{N), resulting in the momentum dependent components 
c“(p,a;o). By an analysis similar to the gluon case, one then gets the boundary 
conditions []64| 


a*c“(p,o) 

c“(p,0) 

c“(p,T) 


0 if /“ G Cn and p = 0, 
0 else, 

0 . 


(4.23) 


The same boundary conditions are valid for the ghost field c. 
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4.1.3 The total action 

Having fixed the gauge, we are now able to write down the action that will be 
used for perturbation theory in the Schrodinger functional. 

The gluonic part of the action may be summarised as 

= ^^w{p)tT{l-U{p)} + Sgi[V,q], (4.24) 

% p 

where the weights w{p) are chosen such that the action is 0(a) improved, i.e. 
w{pt) = Ct and w{ps) = Cg/2 for the time-like and space-like plaquettes at the 
boundaries. 

The ghosts contribute to the action with S'pp given in ( |4.21| ). The fermionic part 
of the total action is given by the Dirac-Wilson action including the volume and 
boundary improvement terms, 

Sf[V, = a^'^ip{x){D + 5D + niQ)ip{x). (4.25) 

X 

Apart from these contributions, there is an additional part of the action arising 
from the change of integration variables. Before, the gauge field integral was 
formulated using the measure D\U]. Now, we want to express the Schrodinger 
functional as an integral over the gluon fields q. For this purpose, it is useful to 
use 

m = n 

as a new measure. It may be obtained from the old measure by 

D[t/] = (4.27) 

where, in the case of SU{3), the measure part of the action is given by 

= f I] + 0(9^). (4.28) 

x,fi,a 

In the case of SU{2), one would have to replace the factor | by A. go the 
additional contribution to the action due to the change of integration variables 
turns out to be of order g^. For the later calculations in this thesis, the total 
action will only be needed up to order go, so may be ignored here. It has, 
however, to be taken into account in the computation of the running coupling at 
2-Ioop order. 

The total action is the sum of all contributions, 

S'tot[F,g, c, c,'0,'0] = SG[V,q] + S'Fp[V,g, c, c] Fp)!/, g,-0,-0] + S^[q]. (4.29) 
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The Schrodinger functional now becomes 




D [g] D [c] D [c] D [ip] D [V^] . 


(4.30) 


Let now (9 be a product of link variables. (We assume that O does not 
depend on the quark helds, which is the case after integrating over them.) The 
expectation value is then to be computed using the probability density (p.61|) . 
Writing the quark determinant as an integral over Grassmann variables, one gets 


{0)g = ^j (4.31) 

As a product of link variables, O may be written as a series in go, 

O = + goO^^^ + + 0(g3), (4.32) 

where (9^"^ contains products of n gluon fields q^. Since is merely a constant, 
the functional integral becomes trivial, 

= O®. (4.33) 

One thus gets the expansion 

{0>o = 0<») + 90 + si + 0(9“). (4.34) 

However, this is not yet the whole expansion, because the total action and thus 
the expectation value of still depends on the bare coupling. In order to get 
the correct series, one has to expand the action too, 

'S'tot = —‘S'tot ^ + “S'® + fi'O'S'tot + 5'o‘S'tot + O(5'o)- (4.35) 

9o 


f _2) 

The term S'tot comes from the gluon action and contains only the background 
field. 


c.(-2) 

*Jtot 


5^tr{l -l/(p)} 


TT^ 

12^[sin2(7) + 2sin2(7/2)], 


(4.36) 


with 


^ LT ^ 3 


(4.37) 


/_2') 

Since S'tot ^ does not depend on the variables integrated over in the functional 
integral, it cancels in all expectation values and may be neglected here. With 
this expansion, the exponential factor in the functional integral becomes 


T - Stot 


1 ~ 5'o‘S'tot + 9o ( 2‘S’tot 


(1)2 




( 2 ) 

tot 


_c(0) 
g ‘^tot 


(4.38) 
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Let now (. • .)o denote the gauge field average taken at zero bare coupling, which 
means that this average is computed using S'® instead of the whole action. 


Inserting ( 4.38 ) into the functional integrals and integrating out the quark fields 
yields 


{0)q = C>® + 


9l 


(c>d) 


c(l) 

^tot 

Jl 


F/ 0_ 


+ 0(5'i 


(4.39) 


Note that no terms proportional to go or occur, because this would involve 
an integral over an odd number of gluon fields, giving zero. Another observation 
worth mentioning is that S'^® only is non-zero due the presence of the background 
field. Later, this term will lead to some Feynman diagrams that do not occur in 
the case of a vanishing background field. 


4.2 Perturbative expansion of the coupling 


4.2.1 The coupling at 1— and 2—loop order 


The Schrbdinger functional coupling g"^ can be expanded in powers of the bare 
coupling go, 

f = gl+ Pi{L/a)g^ + p 2 {L/a)gl + 0(^^). (4.40) 

At tree level, both couplings coincide due to the normalisation of g'^. The 1-loop 
coefficient pi may be decomposed in an Af dependent and an Af independent 
part, 

Pi= Pio+PiiNi. (4.41) 


Both pio and pn are easy to calculate by expanding the total action. From the 
quadratic parts of the gluon, ghost, and quark actions one gets the inverse of 
the propagators. The coefficients pio and pn are obtained by differentiating the 
determinants of the inverse propagators with respect to the parameter p. The 
coefficient pio has been computed in (and in for the case of SU{2)), the 
calculation of pn can be found in . 


At 2-loop level, one has a decomposition up to order Ay 

P2 = P20 + P21 Af + P22Af^ 


(4.42) 


Here, the computation is more complicated and requires the calculation of Feyn¬ 
man diagrams. The coefficient p 2 o for the quenched case has been computed 
in 1^, and in for the case of SU{2). The remaining coefficients p 2 i and P 22 
can be found in [R5 . 


4.2.2 The step scaling function and its lattice artefacts at 
1— and 2—loop order 

Having expanded the coupling, one easily gets the expansion of the step scaling 
function, since it is given by the difference between the coupling at length scale 
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2L and L at the required order of perturbation theory. This means one has to 
calculate the difference 


Apij{L/a) =pij{2L/a) -pij{L/a). 


(4.43) 


Also the deviation 6 of the step scaling function from its continuum limit (^.731) 
can be expanded in perturbation theory, 


S[u, cl/ L') — [(5io + (5iiiVf] u + [(520 T <52iiVf + 522 ^/"^ + 0('u^ 


(4.44) 


The continuum limit of the step scaling function is determined by the beta func¬ 
tion. In order to get the coefficients 6 ij, one has to decompose the coefficients of 
the beta function, (|2.26| ) and (|2.27|) , into the iVf independent parts feoo and 6io 
and the Nf dependent parts 6oi and fen. The coefficients 6 ij are then given by 


<5io 

— ^Pio ~ 2feoo ln(2). 

(4.45) 

fell 

= Apii - 2feoi ln(2). 

(4.46) 

fe20 

= Ap 2 o - 2feio ln(2) - 2Apio(pio + feoo ln(2)). 

(4.47) 

fe21 

= Ap 2 i — 2feii ln(2) — 2Apii(pio -|- feoo lii(2)) 



-2Apio(pii feoi ln(2)). 

(4.48) 

fe22 

= Ap 22 - 2Apii(pii -F feoi ln(2)). 

(4.49) 

loop coefficients fei ^ were first listed in |66| , while the 2-loop coefficients fe 2 i 


The 1 

have been estimated in P4 
two-loop coefficient 62 contains the critical quark mass at 


However, at vanishing renormalised quark mass, the 

1-loop order, which 
In order to do a 


only in the continuum limit reaches the value used in [34 


precise calculation, one has to specify the zero mass condition with the cutoff in 


place. So the results of [^ could only give a first idea of the size of the cutoff 
effects. For more accurate results, the critical quark mass has to be expanded up 
to 1-loop order at finite a/L and inserted in the formulae for 52j in 0] 
will happen in chapter 


This 


According to the estimation in [34], the cutoff effects seem to be small. In or¬ 


der to compute the step scaling function non-perturbatively, one has to simulate 
a sequence of lattice pairs with decreasing lattice spacing and fixed coupling u 
and extrapolate the Monte Carlo data to the continuum limit. In this procedure, 
the perturbative expansion of (5 (m, a/L) may be used to remove the cutoff effects 
up to 2-loop order from the non-perturbative values of the step scaling function 
S(M,a/L) 1^. 





























Chapter 5 

The current quark mass 


5.1 The PC AC relation 

Several ways to define a renormalised quark mass have been sketched in chapter ^ 
For the Schrodinger functional, a suitable dehnition has been given in |Q, using 
the partial conservation of the isovector axial current. The axial current and 
density are dehned by 


= V'(a^)7/.75^rXx), 

(5.1) 

P“(a:) = ?/5(a:)75^rXa;). 

(5.2) 


As already mentioned in chapter a dehcit of Wilson fermions is that the chiral 
symmetry of the theory is explicitly broken and only restored in the continuum 
limit. Thus only in the continuum limit, the isovector axial current A“(a:) will 
satisfy the PCAC relation 


d^Alix) = 2mP“(x), 


(5.3) 


while at hnite lattice spacing, this relation will be violated by terms of order a. 

In order to make the axial current converge faster towards its continuum limit, 
one may apply Symanzik’s improvement programme by adding an appropriate 
improvement term 

(Ai)“(a;) = A“(x) + aM“(x). (5.4) 

The improvement term turns out to be 


M“(x) = CA-(a, + a;)P“(x). 


(5.5) 


This 0(a) correction is proportional to the improvement coefficient ca, which is 

ca{9o) = -0.00756(1)(?2 + (5.6) 
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to 1-loop order of perturbation theory ||^. The axial density P“ can be shown 
to converge to its continuum limit with a rate proportional to a^, hence it does 
not need to be improved. 

In the case of a non-vanishing physical quark mass, the quantities defined 
this way are still not fully improved. In order to improve ^4“ and completely, 
one still has to subtract a mass dependent counter-term. However, this amounts 
to a mass dependent multiplicative renormalisation, so it seems more natural to 
include this factor in the dehnition of the renormalised quantities. 

In general, the quark mass will get an additive renormalisation, which means 
that, in the plane of bare parameters, there will be a critical line 


mo = m.e(go), 


( 5 . 7 ) 


where the renormalised quark mass vanishes. For convenience, one may define 
the subtracted mass 

fnq = — fnc- (5.8) 

Then this mass will only have to be renormalised multiplicatively. 

Including the already mentioned mass dependent factors needed for improvement, 
the renormalised axial current and density may now be written as 


— ^a( 1 + acA-(9^ + 9*)P“}, 

(Pr)“ = Zp(l + 6pamq)P“. 


( 5 . 9 ) 

( 5 . 10 ) 


Here, &a and bp cancel mass dependent cutoff effects. At tree level one has ^ 


5 ( 0 ) _ do) _ . 


( 5 . 11 ) 


Since all calculations in this thesis are done at vanishing renormalised mass, which 
is equivalent to setting ruq = 0, these coefficients are not needed here. In the 
following, they are therefore ignored. 

The renormalisation constant Za would be equal to one in the continuum. 
This is a consequence of the SU{Nf)p x SU{Nf)p symmetry, leading to chiral 
ward identities, which may be used to normalise the currents |]6^ , [69[| . In the 
regularised theory, this symmetry is violated by terms of order a. This means 
that the PCAC relation and the chiral ward identities are only valid up to 0(a) 
corrections, resulting in a renormalisation of the axial current. It stays, however, 
finite and scale independent. Up to 1-loop order of perturbation theory, one 
gets 1^, ItTI 

Za(^/o) = 1 + 4'^1?o + 0(^?o'), (5.12) 


with 


Zy = -0.087344(1) X Cp, 


(5,13) 
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In contrast to the renornialisation constant Zp is scale dependent and thus 
responsible for the running of the renormalised mass. This issue will be addressed 
in chapter 

A renormalised quark mass mp may now be dehned as the proportionality 
constant in 

+ d,){A^)l(x)o\ = 2mR ((P^mo) + O(o^) (5.14) 

for any product O of renormalised improved helds located at a non-zero distance 
from each other and from x. The lattice artefacts of order depend on the 
choice of O. 


5.2 The current mass and its lattice artefacts 


In order to compute the renormalised mass, one needs to choose an operator O. 
Using this operator, one then has to construct bare correlation functions contain¬ 
ing the axial current and density and then renormalise both these quantities and 
the operator O. 

One possible choice is to introduce the bare correlation functions 


IaM 


fp(xo) 


J^UAS(xK(y)7^^T'‘c(z)\, 


y,z 


-a® X ( P“(a;)C(y)75xr“C(z) ) , 


r“C(z)^, 

(5.16) 

r“c(z)y 

(5.16) 


where C (x) is the functional derivative with respect to the boundary quark helds 
at Xo = 0 dehned in ( p.27|) and (|3.28|) . These correlation functions are propor¬ 
tional to the probability amplitude that a quark antiquark pair created at Xq = 0 
propagates into the interior of the lattice and annihilates at the point x. This 
situation may be depicted by hgure |0| . 

Now one can dehne the xq dependent current quark mass 


m(xo) 


^((9o -F (9 o)/a(xo) -t- CAagogo/p(xo) 
2/p(xo) 


(5.17) 


which, taken at a certain xq, may serve as an unrenormalised quark mass. A 
convenient choice is to take m in the centre of the lattice, i.e. 


J m for even T/a, 

I 7 for odd T/a. 


(5.18) 


This current quark mass will, of course, depend on the lattice size L/a. It is 
then straightforward to show that the renormalised quark mass dehned as the 








40 


The current quark mass 


Xo= 


Xo= 



LxLxL 

Figure 5.1: The correlation functions /a and /p 

proportionality factor in the PCAC relation is given by 

Za(1 + bAamq) 


rriK = 


-mi + 0(a). 


(5.19) 


Zp{l + b^am^ 

One may of course take different operators O which will give quark masses dif¬ 
fering from mi by terms that are expected to be of order a^. These 0(a^)-terms 
may serve as a direct check of the size of the lattice effects. 

A possible choice is to take the boundary quark helds at Xq = T instead of 
Xq = 0, i.e. to use the correlation functions 


/a(T-xo) = /A“(x)C'(y)75^r“C'(z)y 

/p(T-xo) = /p“(a;)C'(y)75^r“C'(z) 


y.z 


(5.20) 

(5.21) 


instead of /a and /p. Let m! be dehned in the same way as mi, but with and 
/p instead of /a and /p. Now, a new unrenormalised mass may be dehned by 


m2 = + "^0- 


Then the lattice artefact 


d{L/a) = m^^L/a) — mi{L/a) 


(5.22) 


(5.23) 


should be of order a^. 

A different possibility to estimate the size of the lattice effects is to take the 
difference 

e{L / a) = mi{2L / a) — mi{L / a). ( 5 - 24 ) 
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One aim of this work is to compute these lattice artefacts up to 1-loop order 
of perturbation theory. To this end, the 1-loop expansion of the correlation 
functions /a and /p is needed. 


5.3 1—loop expansion of f\ and /p 


For convenience, the calculation in this section is done in lattice units, i.e. a = 1. 

The principle ideas of the expansion of /a and /p have been outlined in 
a more detailed discussion has been given in [^] (and in ||^ for the case of a 
vanishing background held). In this section, the details of the calculation shall 


be presented following the lines of . 


5.3.1 Preliminaries 

In order to expand f\ and /p, they have to be brought into a form more suitable 
for perturbation theory. Applying Wick’s theorem on (|5.15|) and (|5.16|), one gets 


/a(xo) = 




fpi^o) = 


z)V'(x)]F7o75[V’(a^ 

z)i>{x)]F'y5[i>{x 




where the trace is to be taken over Dirac and colour indices only. 
Let 'ijjci be a solution of the Dirac equation, 

{D + 6D + mo)'ijjci{x) = 0. 

Then it can be shown that 

Helix) 


(a;)C(y)]F = 


So using the matrix 


H{X) = 


Hiy) ' 

Hcljx) 

<5p(y) ’ 


/a and /p may be written as 


fAixo) = --{Tt{H{x)^'JoH{x)})^, 

/p(i„) = 1 ■ 

The matrix H has an expansion 

Hix) = H^^\x) + g,H^^\x) + glH^^\x) + 0(^7^), 


(5.25) 

(5.26) 


(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 


(5.32) 
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resulting in the 1 -loop expansion of /a, 
1 


x^'yaH^^^x) 


Tt 
- /Xr 


c(l) 

*Jtot 




X, 




( 1 ) 

tot 


F/ 0 


(5.33) 


and analogously for fp with —1 instead of 70 . So, in order to compute the 1-loop 
coefficients of /a and fp, one has to expand Stot up to order qq and 


H{x) = Ct'^S{x,y)U{y -6,0) 


-ip_ 


+ lyo=i 


(5.34) 


up to order g^. Here, S{x, y) denotes the quark propagator, which is the inverse of 
the Dirac-Wilson operator {D + 5D + mo). So one has to expand the propagator 
S{x,y) and the link variables. Since the overall factor ct cancels in the dehnition 
of the current mass, it may be ignored here. There will however be a contribution 
proportional to due to the boundary counter-term in the quark action, leading 
to a corresponding term in the gQ coefficient in the expansion of the propagator. 

Following the lines of |M, the calculation is done in momentum space. First, 


the gluon helds have to be decomposed in a basis of the Lie-algebra of SU (3) 


= '^q'^{x)P, (5.35) 

a 

with the basis J“ chosen such that the star operation dehned in appendix ^ and 
the covariant derivative act diagonally: 


cosh Gofc ^ “ 
sinh Gok ^ 
{Dkf){x) 

{Dlf){x) 


CJ\ 

SJ\ 

^ [fla(xo)/“(x + k)- P{x) 


P, 




VLa{Xo)*r{x 



(5.36) 

(5.37) 

(5.38) 

(5.39) 


The /“ are explicitly given in appendix The coefficients f2a(a:o) can be 
parametrised as 

f^a(xo) = (5.40) 

with the given in appendix]^. Then one may introduce the Fourier trans¬ 

formed helds 


x) 


P 


(5.41) 
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x) 


1 

Z3 





P 

Furthermore, the 0^ may be used to define the quantities 


4 (p,xo) = 2 sin 

^{pk + M^o)^ 

Cfc(p,Xo) = 2 cos 

^{pk + 4>a{Xo)^ 


and 


(5.42) 


(5.43) 

(5.44) 


fla = (C„ - (5.45) 

which will be used to make the expressions for the propagators and vertices more 
compact. Here, Ca and Sa are the components of cosh Gok and sinh Gofc according 
to ( |5.36| ) and (|5.37D . They and the Ra are explicitly given in appendix ^ 

The same is done for the ghost fields 

P pa 

P pa 

and also the quark fields are Fourier transformed 

0(x) = ^^e*P^0(p,Xo), (5.48) 

p 

ii{x) = ^^e*P^0(p,xo). (5.49) 

p 


( 5 . 46 ) 

( 5 . 47 ) 


5.3.2 The propagators 

In momentum space, the quadratic part of the pure gluonic action takes the form 


^g^ = I^Y Y Y€^~P^^o)K^uiP',Xo,yo)q^{p,yo), ( 5 . 50 ) 


T-l 


2L3 


P a:o,yo=0 a 


where 


Kki{p-,xo,yo) = 5^o,yo 4/s“(p,xo) - Sfc(p,xo)s“(p,(ro)(l - Aq) 


R^kl 


‘^CaSxQ^yQ Rai^^xo+l,yo T —l,yo) 


(5.51) 
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^fco(p; ^ 0 , yo) 


^^ofc(p;^o,i/o) 
^oo(p; ^ 0 , 2 / 0 ) 


%Rn 


<5xo,yo4(P,^0 + 1) - <5x0-1,yo4(P, 2 / 0 ) 


-iAo4(P,31o) 


<5x0-1,: 


■'Xo,J/o "xo-i,j/o 

-Ra*5a;o,yo s“(p,a;o) ■ s“(p,xo + 1) 

+Ao 25x0,1/0 ~ <5x0+1, 1/0 ~ <5x0-1,1/0 
Ao5xo,yo <5x0, 0(1 Xa<5p,o) ^xq,T —1 


(5.52) 

(5.53) 


(5.54) 


where Xa = 1 for a = 3, 8 and Xa = 0 otherwise. Then the free gluon propagator 
is given by 

(?p(p,a;o)gv(p'>2/o))o = <5fea^^<5p+p'T>“^(p;xo,i/o), (5.55) 

where iA“ is the inverse of K°^. 

The ghost action takes the form 

CXD ^ 

n=0 

where the quadratic part is given by 

^o)^“(p; So, to)c“(p, to)- (5.57) 

p So,to “ 


For a 7 ^ 3, 8 one has 


F“(p; So, to) = Sso,to [2 + s“(p, so)^] - Sso+i,to - <5so-i,to- (5-58) 

The computation aimed at here only involves closed loops of charged ghost helds 
coupled to a neutral gluon. Hence the quadratic ghost action for a = 3, 8 is not 
needed here. Now we get the free ghost propagator 

(c“(p, So)b^(p', to))o = 56aF^5p+p'T>“(p; sq, to), (5.59) 

where iA“ is the inverse of F“. 

The free part of the quark action is given by 

4”' = ^EE i5(-p, xo)D{p] xo, yo)'i/j{p, yo), (5.60) 

P xo,yo 

where D is the improved Dirac-Wilson operator [D + 6D + mo) at lowest order 
of go, namely 


D{p,Xo,yo) -P—<5xo+l,i/o 4"-®(P ,^ 0 )^x 0 , 1/0 -^+<5x0-l,yo’ 


( 5 . 61 ) 
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with P'1 = pk + 9/L and 

B{p,xo) = 4 + mo-^ 

k *- 

+ihr7o^7fc. 


i(l + 7k)e-^^’^Vixoy + ^(1 - 7;^)e*^’'=f^(a;o) 


(5.62) 


Here hf is a diagonal matrix in colour space with the elements 


Haa = -xciw sin£’„. 


(5.63) 


The quark propagator S is now defined as the inverse of the operator D. 


5.3.3 The vertices 

The expansion of the pure gluonic action gives the triple gluon vertex. At 1-loop 
order we have 

^ ^p(qi + q 2 + qs) 

qi.q2,q3 

■ X] X] X] i^/ri)?2/^‘|(Hnq2,q3;H,t2,t3) 

M1)M2>A‘3 tl,t2,t3 0.1,0.2, 0.3 

(5-64) 

3 

where (5p denotes the periodic delta function, i.e. the delta function modulo 27r. 
In this expression, a term proportional to is missing, which will be treated 
later. In order to make the somewhat complicated expressions for the vertex 
V slightly more compact, we introduce the following notations for traces and 
permutations of the basis of the Lie algebra and the background field: 

Cabc = -2i trJ“[J^T], 

eabc = -2i ri^F - 6-*^/"/'’/“). 

Furthermore, we use the shorthand notations 

— ^^{Pk+Qk+rk) 

and 0), = docj)c{xo). Since all non-zero (j)c are linear in xq, cj)'^ does not depend on 

Xq. 


(5.65) 

(5.66) 

(5.67) 
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Using these notations, one gets for the vertex 

^Zm(P) ^ 0 ) ^o) Cabc^so,to^so,uo^^klEkEm 

■ sin^lpm - 4>a{so) -qm + (j)b{so)]cl{-r, Sq) + 2 permutationsj 
% % ( 

^ 2 ~l~ 0 fo(^o) “t“ 0 c(^o)]^A;/m-^fc “t“ ^bac)^so+l5^iO 

-{ecab + ec6a)(5so-i,no] + 2 pemiutations|. (5.68) 

Here, “2 permutations” is meant such that one has to take cyclic permutations of 
k,l,m and a, b, c, respectively. The vertex is thus made totally symmetric under 
interchange of labels. Then the other parts of the vertex become 

V))oo''(p, q, r; So, to, Uo) = -^(55o,to'^so,«o ^{((^abc + 

k 

— (scab + + 2 permutationsj, (5.69) 


^fc“o‘'(p, q, r; So, to, ^o) = -^bkiEk 
1 


-K,toK,uo{{eabc + - {e,ab + 


r*^soiio*^.50—i,uo(^cafe “h ecbQi)e ^'^'^s^( r, mq) 




e 2 n 


^cab^ 




So,Uq 


+ 1 permutationj), 


(5.70) 


2 

^oofc'(P) q> r; So, to, Uo) = --Sso,to 

■{-lK,uo{{ecab + e,).,)e-5h.-^d«o)] + } 

— “t“ ^6ac)Cfc( '^o) 

— p ^ J 

[^acb'-^so,UQ ^abc'-^SQ,UQ — l\ 

+ 1 permutationj j. (5-71) 

The gluon-ghost vertex is given by the 1-loop expansion of the ghost action. 
Here one gets 

‘^gh = <5p(p + p' + q) 

p,p',q so,to,uo /7 

^(-p', So)7^“'’''(p', P, q; So, to, Uo)c\-p, to)q))(-q, Mo)- (5.72) 

a,b,c 
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As already mentioned, only the vertex of a neutral gluon (i.e. c = 3, 8) and 
charged ghosts is needed here. In this case, the vertex becomes 


(5,73) 


^0 (p ) P) Q)'SQ; ’^o) [^so+l,io A (5so,io] 

—l,no [^SQiio A li^o] M 


''(p',P,q; So,^0,Mo) = Cabc5so,toSso,uoCk{-P,So) 

■ sin-[-Pk -qk + M^o) + 0 c(so)]- (5.74) 

The quark action has to be expanded up to order gQ to get the quark-quark- 
gluon and the 2 quark - 2 gluon vertices. Both vertices have two terms, one 
coming from the Wilson part and one coming from the Sheikholeslami-Wohlert 
part of the action. 

The Wilson part of the action has the expansion 


C _ c(0) I \ ^ „n cjn) 

>JF,Wilson — Wilson i/O ^F,Wilson’ 


(n) 


(5.75) 


n=l 


where 4°wiison 


has been given in (|5.6CI|) and the n-th order expansion coefficient 
may be written as 


q{n) 

‘^F,Wilson 


1 

l3 


n +1 




p,p',qi,...,q„ 


i=i 


fl Xo,yo,ZO 0,lr--:0.n 


^/’(-p', aio)V;r''"“"(s; xo, Vo, zo)ij{-p, yo) JJ (-q^-, zq), (5.76) 

i=i 


with s = 2 (P^ ~ P)- Using the notation 

jai...an _ ^ \ ^ J'7(q.i) jcr{an) 

n\ ^ ’ 

perms a 

the time components of the vertex V become 


(5.77) 


Vo^^-^-(s;xo,yo,zo) = 

+ {-l)'^P+Sxo-l,yoSyo,zo^. (5.78) 

In contrast to the time components, the space components 

V,^^-^-is;xo,yo,zo) = -/“^•••“’^^4o.,o<5.o,.o{u^r-“"(s,a;o)(l-7it) 

-1/ 


+ (-l)"W'“--““(s,a:„)-‘(l+ 7 t) 


(5.79) 
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depend on the background field via 


ITf •••“"(s,xo) = l/(xo)exp(^z[s+ + ^ ^^o)]) • (5.80) 


While the Sheikholeslami-Wohlert term (|2.14|) does not contribute to the free 
action (^.OOp, it has to be taken into account for higher orders, 


= Csw(<?o) V (5.81) 

n\ 

n=l 

where the n-th order expansion coefficient is given by 

= (^)”" E ^p(p+p' + E<i,) E E E 

^(-p', a;o)*5p(:::p:(qi, . . . , q^; ^10, . . . , ^no)V'(-P, Xo) 

n 

(5-82) 

i=i 

For n = 1, the vertex S may easily be written down as 


S'o(q;a;o,^o) = ^ ^oic“(-q, Xq) 

j 

■ ^(*^xo,zo “t“ —l, 2 o) ,COs£^]-Sj( q, Xq) 

-(5x0,20 -'^xo-i,2o)[^“>sin^]c“(-q,a;o)), (5.83) 


Ski(llXo,Zo) = ^aofc(^(5xo+i,20 - 5x0-1,2o){^“>cos^}Cfc(-q,zo) 

(5x0+1,20 ‘^^xo,zo 5x0-1, zo) [.^ q^ ^o)^ 

-^/“5xo,zoCfc(-q, ^o) (“^’ ^o)c“(-q, ^o). (5.84) 

j 

For n = 2, the vertex is much more complicated. Fortunately, it is not needed 
completely for the calculation of /a and /p at one loop order. The only thing 
which really is needed is at order gQ, which is given by 


{^jk{x)) = 0 {g^), 


( 5 . 85 ) 
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and 

q c ^ 

2i cos[£: -^{qk + 4>c{xo + 1))] (^oo(q; ^ 0 ,2:o)Sfc(q) ^0 + 1) 

+2*£’fco(q;^o + l,a;o)) 

-2i cos[£: + ^(gfc + Mxq - 1))] ^0 -^,xo- l)4(q, ^0 - 1) 

-2*^fco(q; ^0 - l,a;o - 1)) 

+i sin £ (^Dlj^{q; Xq + 1, Xq + 1) + 2D^fc(q; Xq, Xq) 

+^fcfc(q; xo -i,xo- 1)) 

-2i sin[£ - ^0^] Xq + 1, Xq) + 4^fcfc(q; Xq - 1, Xq) 

+i4(q, xo + i)i^fco(q; ^0, 2 : 0 ) - *Sfc(q, - i)^fco(q; ^o, xo - 1)) | 
+O(^7o'). (5.86) 

The complete vertices may now be computed by combining the Wilson and 
Sheikholeslami-Wohlert parts. The quark-quark-gluon vertex is then given by 

V'^“(p',P,q;so,^o,Mo) = V^{s; So,to,Uo) + So,Uo)6,,^to, (5-87) 

and the 2 quark - 2 gluon vertex by 

V),t(p',P,q,q';5o,to,Mo,Mo) = V^\s]So,to,Uo)Sf,Juo,u'o 

+c^w(q> q'; ^o, ^o, UQ)^so,to- (5.88) 


5.3.4 The diagrams 


Having calculated the propagators and vertices, one may now expand /a and /p. 
The calculation can be made slightly more general by dehning the function 

f{r-,p,x„) = ^e‘»'>'-">i(Tr{F+rp_t;(s-o,o)S(z,i) 
y>z 


■r^(x,2/)t/(i/-0,0)-^} 

According to (p.30| ), (|OT|) , and (|5.34| ), one has 


yo=zo=^ 


(5.89) 


/a(xo) = qV(7o75;0,Xo), 

/p(Xo) = CtV(75;0,Xo). 


(5.90) 

(5.91) 
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Xo=0 • 


Figure 5.2: Diagram for /a and fp at tree level. The dotted lines denote the link 
variables between xq = 0 and xq = a, and the cross represents the insertion of 
the axial current or density. 


Calculating the function / at tree level only amounts to taking the tree level 
values of the quark propagators and Fourier transforming everything, leading to 

/(F;p,Xo)(°) = ^Tr|p+FP_F(p;l,a;o)FF(p;xo,l)}. (5.92) 

Since the quark propagator is diagonal in colour space, f(T;p,xoY^'> be 
decomposed into colour components 

N 

/(F;p,Xo)^°^ = 5Z/(r;p,Xo)i°\ (5.93) 

a=l 


with the colour components given by 

/(r;p,Xo)L°^ = ^tr|p+FP_F(p;l,Xo)ar^(p;xo,l)a}, (5.94) 

where tr denotes the trace over the Dirac indices only. Thus one gets a tree level 
coefficient which can be depicted by the diagram in hgure |5.2| . 

For the 1-loop expansion, the propagators and link variables in (|5.89| ) have 
to be expanded according to 

S{x,y) = S^°\x,y) + goS^^\x,y) + glS^^\x,y) + O(^o)> (5-95) 

and 

U{x,ii) = exp{goqf,{x)}V{x,iJ,) 

= +9o(l^,{x)+ ^gQqf,{xY+ 0{g^)"^V{x,y). (5.96) 
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In order to compute / at 1-loop order, one has to insert these expansions, gather 
all contributions of order and contract the gluon helds to gluon propagators. 
Apart from these contributions, there are improvement terms proportional to the 
1-loop improvement coefficients. In order to compute / at mi = 0, one has to 
do the calculation at the critical quark mass me, which has an expansion 

me = m^e"’ + + O(5'o)- (5.97) 

The 1-loop coefficient of / then gets a contribution proportional to m^c^. Thus 
we get the sum 


/(r;p,Xo)^^^ = + ci^f{T;p,Xo) 


( 1 ) 

V 


+cl^V(r;p,xo)Fb + ci"7(r;p,a;o) 


(1) , Ji). 


Gb 


+m, 


( 1 ). 


d 


dmi 


■/(r;p,a;o) 


(0) 


(5.98) 


Here, all expansion coefficients have to be calculated at mo = m^\ which will be 
obtained by a tree level calculation explained in section |^ . For /a, one gets the 
additional contribution proportional to 

The contributions /(T; p, a;o)n ^ may be depicted by the diagrams in hgure | 
The second order contributions of the link variables give diagrams la and lb. 


/(riP,a:„)<L> = -~^/(rip,a:„)f^ J]B5,(q;0,0)C 

a q a 


(5.99) 


where 

v-' )aa * 

Diagram lb is simply given by 


f(^',P,xo)u = /(r;p,xo)S- 


(5.100) 

(5.101) 


In the following, we dehne r = p + q for convenience. Then the contraction of 
the hrst order terms of both link variables yields diagram 2, 

/(F; p, = T ^ ^ /(F; r, Z)S„(q; 0,0)C;. (5.102) 

Q; q a 

Obviously, there are four different ways to combine the hrst order terms of the 
link variables with the hrst order terms of the propagators. Contracting the hrst 
order term of the link variables with the hrst order terms of the propagators on 
the same side gives diagrams 3a and 3b, 

/(r;p,x„)<i> = E E^S„(q;0.«o)Tv{rp+rp_s(pii,x„) 

q n SQ,to,uo a 

■rS'(p;xo,So)ld“(p, -r,q; So,to,Uo)S{r-,to, 1)|, (5.103) 
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Figure 5.3: Diagrams contributing to /a( 2 : 0 ) aud fp{xo) at 1-loop order of per¬ 
turbation theory both at vanishing and non-vanishing background Geld. 
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/(r;p,i„)« = E E®^o(‘l!“»>°)Tr{r-P 4 -rP_S(r;l,s„) 

q so,to,uo a 

-p, -q;so,to,Uo)S{p;to,Xo)TS{p;xo, 1)|, (5.104) 

while combining them with the hrst order terms of the propagators on the oppo¬ 
site side results in diagrams 4a and 4b, 

/(r;p,io)S = ^;^EE E E^S 4 'i'°''“«)^'{^“-P+r^’-'S(pa.so) 

q fi so,to,uo a 

'U“(p, -r, q; So, to, Uo)S(r; to, a:o)rS(r; ij, 1)|, (5.105) 

q /I so,to,uo a 

■rS(r;xo,So)V^^(r, -p, -q;so,to,Uo)S(p;to, 1)|. (5.106) 

The second order terms of the quark propagators consist of two parts: One 
containing two quark-quark-gluon vertices with a quark line in between, and 
one containing a 2 quark-2 gluon vertex. The part with the two vertices leads to 
diagrams 5a and 5b, 

f(r;p,f:o)ii> = ~EE E E 

q /J.,iy so,to,uo sl^,t'g,u'g a 

■S'(p; 1, a;o)rS'(p; Xq, So)14“(p, -r, q; Sq, to, Mo)<S'(r; to, 4) 

■'C(r, -p, -q; 4, ^ 0 , “o)>^(p; ^o, i)}, (5.107) 

/(riP,s:o)S’ = ^;^EE E E E “i. «o)Tr{ AFP- 

q 11,V so,to,uo SQ,tg,Ug a 

■S'(p; 1, So)14“(p, -r, q; So, to, Mo)5'(r; to, 4) 

-p, -q; 4,4, 4)>^(p; 4,2 ;o)r,S(p; xq, i)|, ( 5 . 108 ) 

while the part containing the 2 quark-2 gluon vertex gives diagrams 6a and 6b, 

/(riP,s:o)iL’ = -~EEEEE^“s(q!“i.“»)^4+rp- 

q fi,i/ so,touo,UQ a 

■S(p; l,Xo)rS(p;xo,So)Vp(p, -p, q, -q; So,to,Mo,Mo) 
■^(p;to,l)}, (5.109) 
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/(r; P. <■ “»)Tr{p+rF_ 

q fl^U SQ^tQ UO,Uq Q. 

■S'(p; 1, So)V'“(p, -p, q, -q; Sq, to, ^o, Wo)5'(p; to, a^o) 
•rS'(p;xo, 1)|. (5.110) 

In the improved theory, these diagrams do not give the total contribution of the 
second order terms of the quark propagators. There are two more parts, one 
proportional to and one proportional to Csw^ These contributions will be 
dealt with separately. Finally, the combination of the hrst order terms of both 
propagators results in diagram 7, 

/(r;p,x„)t> = i;^EE E E 

q 11,u so,to,uo s'Q,t'Q,u'Q a 

■5'(p; 1, So)V)f(p, -r, q; Sq, to, Uo)S{r; to, a;o) 

■r^(r; xo, so)iC(r, -p, -q; 4,4,4)>^(p; 4, !)}• (5.111) 


So far, the calculation has been completely analogous to the case of a vanishing 
background field. However, with the non-zero background field considered here, 
one has to take the terms in (|5.33|) containing S^l into account. This means one 
has to contract the first order terms of the link variables and the propagators with 
the hrst order terms of the total action, leading to the diagrams in hgure 
Starting with the hrst order term of a propagator, we get 


/(r;p,io)S = -lEEEE^J,.(Oi“i'“»)^{^+r^- 

11 , 1 / So,to uo,Ug a, 

■S{p; 1, a;o)rS'(p; xq, so)lC(p, -p, 0; sq, to, 4) 
■^(p;to,l)}T“(«o), (5.112) 

where T“ denotes the sum of the closed gluon, ghost, and quark loops, 

4^(“o) = 4^,gluon(“o) + 4^,ghost(“o) + 4^, quark (“o)- (5.113) 


Note that the time component Tq vanishes due to CP invariance. The gluon loop 
is given by 

1 1 X - 

4^“gluonK) = 0; So, to, Mo)^4(q; to. So), 

q u,p so,to c 

(5.114) 

the ghost loop by 

ghost (“o) = Yi EEE 0; So, to, Mo)^ 4 q; ^O, so), (5.115) 

q so,to C 
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Xo=0 • 

9a 



Xo=0 • 

9b 


- 





Figure 5.4: Tadpole diagrams contributing to fA(xo) and /p(xo) at 1-loop or¬ 
der of perturbation theory. These terms are only present with a non-vanishing 
background Geld. 
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and the quark loop by 

quark (wo) ^3 EE Tr|l/“(q, -q, 0; Sq, to, Mo)5'(q; to, So)}- (5.116) 

q So,to 

Analogously, one gets diagram 8b, 

lj,,u So,to uo,Uq a, 

■A(p; 1, So)K'(P, -P, 0; So, to, Mo)*^(p; to, 2 : 0 ) 
■r^(p;xo,l)}T;(no). (5.117) 

The combination of and the hrst order terms of the link variables yields 
diagrams 9a and 9b, which are given by 

/(rip,i„)S = -l5^5^5^DyoiO,«„)Ti{/“P+rp_s(pii,x„) 

fi iLQ a 

■r^(p;xo,l)}T“(no) (5.118) 

and 

/(r;p,Xo)9l^ = -/(r;p,Xo)£^ (5.119) 

Since they are of opposite sign, they cancel in the sum and may be ignored. 
This leaves diagrams 8a and 8b as the only diagrams depending on the number 
of flavours. Note that these terms are not present in the case of a vanishing 
background held. This means that with a non-vanishing background held, the 
critical quark mass is dependent on iVf at 1-loop order, while with a vanishing 
background held it is not. 

5.3.5 The improvement terms 

It was already mentioned that diagrams 6a and 6b get additional contributions 
from the volume and boundary counter-terms in the quark action. 

The volume term gives corrections proportional to ciw ,which have already been 
taken into account in the vertices, and corrections proportional to ciw which 
still have to be computed. These corrections to the propagator at order Qq are 
imposed by the Sheikholeslami-Wohlert part of the Dirac operator, which at this 
order is 

= -^4w7o 

k 


(5.120) 
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In diagram 6a, this leads to the correction 


T-l 


/(r;p,xo) 


(h = lp(i) 

6 a, V ^ "^sw 


■7o 


EE tr|P+rP_S'(p; l,a;o)arS'(p;a;o,So)a 

a 50=1 

^7fcS'(p;a;o, l)a [ sin£:„. (5.121) 


Using the symmetry of the quark propagator 

5'(p;xo,2/o) = 755'(p;2/o,a;o)W5, 


(5.122) 


one can show that the contribution to diagram 6b is the same. Thus, the total 
contribution proportional to ciw becomes 

/(r;p,Xo)?^ = 2/(r;p,Xo)gv (5-123) 

Another way to get this improvement term is by Taylor expansion. 


/(r;p,a;o)v^ 


d 

dc 

L/L^SW 




(0) 



(5.124) 


So, the volume term could, in principle, also be obtained from the tree level result 
by numerical differentiation. While this procedure turned out not to be suitable 
for the calculation, since it gives the result with less accuracy, it was successfully 
used to check the result obtained according to ( |5.123|) . 

Apart from the volume term, diagrams 6a and 6b get a contribution from the 
boundary counter-term of the quark action, which acts on the propagator as 

6D^^^S{x,y) = c[^'' + 6^^^T-iS{x,y)} . (5.125) 


In diagram 6a, this leads to the contribution 


/(r;p,a;o)6a^b = -^ci"^Tr|P+rP_^(p; 1, a;o)r ^^(p; Xq, l)^(p; 1, 

+^(p;xo,T-l)^(p;T-l,l) 




1) 

(5.126) 


Here again, the corresponding contribution to diagram 6b can be shown to be 
the same. One thus gets 

/(r;p,Xo)i\l = 2/(r;p,Xo)^t^,b- (5.127) 


Strictly speaking, this does not give the whole term of /a and /p proportional to 
Ct, because there is still the overall factor of in (^.90|) and (|5.91|) to be taken 
into account. As stated before, this factor cancels in the current quark mass and 
may be neglected in this calculation. It may, however, be used as a consistency 
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Im1a 


Im1b 


Figure 5.5: Diagrams for f{T ; p, xq) 


( 1 ) 
V ; 


f{r]p,Xo)pl, ^/(r;p,a:o)^°) 


check, since, as an improvement term, (/a/p pb + ^/a7p) approach 


f(0) 


zero m 


the continuum limit, which /T/ppb /a/p separately. 

Both improvement terms calculated so far come from the insertion of 
into the quark line and may thus be depicted by the diagrams in hgure 5.5 . 

Another term that may be depicted by the same diagrams is the derivative 
of / at tree level with respect to the bare mass. To this end, one has to take the 
derivative 

d 

^(p;a;o,|/o) = - ;xo,So)S{p;so,yo) (5.128) 

^ so 


of the quark propagators in (|5.92|) . The two contributions coming from differen¬ 
tiating the two propagators can be shown to be the same, leading to 


^^/(F; p, = - ^ Tr {P+TP.S{p; 1, Xo)F^(p; xq, So)S{p; Sq, 1)} . 

(5.129) 

With a non-vanishing background held, there is also an improvement term 
proportional to c 77 giving an additional contribution to diagram 8 a and 8 b. 
This term arises from a contribution proportional to in the total action, 
contributing to /a and /p via (|5.33|) . Apart from the gluon tadpole, the gluon 
part ( [4.241) of the total action yields the boundary counter-term 

-^Sb = ^ + sin7], (5.130) 


where the parameter 7 is given by 



(5.131) 
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Im2a 



• • 

Im2b 


=0 


T-a 


Xo=0 



ImSa 



.T-a 


Xo=0 


Im3b 


Figure 5.6: Diagrams contributing to the improvement term proportional to 
for /a and /p at 1-loop order 


Contracting this part of the total action with the first order term of the quark 
line gives the contribution 


/(r;p,a;o) 


( 1 ) 

Im2a/Im3a 



E EE Tr|p+rP_S'(p; l,a;o)rS'(p;a;o,so) 

so,to,uo M k 


(p, -p, 0; So, to, Mo)<^(p; to, 1)| mq, 1) 

-1) [sin( 27 ) + sin 7 ], (5.132) 


which may be depicted by diagrams Im2a and ImSa in hgure |5.6| . The corre¬ 
sponding correction depicted by diagrams Im2b and Im3b can be shown to be 
the same, such that the total contribution proportional to becomes 


f{T-p,Xo)^^l = 2/(r;P,Xo)i^2a/Im3a- (5-133) 

Strictly speaking, there are also contributions proportional to 4^^ from the con¬ 
tractions of the time-like link variables at the boundary with the gluons in (5.1301), 
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Im4a Im4b 




Im5a 


Im5b 


Figure 5.7: Further diagrams contributing to the improvement term proportional 
to for /a and /p at 1-loop order 
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depicted by diagrams Im4a, Im4b, Im5a, and Im5b in figure |5.7| . But, like the 
diagrams 9a and 9b, these contributions are of opposite sign and thus cancel in 
the sum. 

The c[^^-term was neglected in the computation of ciw in |^. This is permis¬ 
sible, because in the current mass, the boundary counter-terms only contribute 
at order a?. This is due to the fact that (^.14 ) remains valid even if the boundary 
counter-terms 6SG,h and 55'F,b are dropped. Further explanations concerning this 

Furthermore, this term was missing in the perturba- 
To which extent this term changes 


matter can be found in 54 


five results for the discretisation errors in 27 


the results published there will be discussed in chapter |^. 

A further improvement term only present in /a arises from the 0(a) correction 
of the axial current (|5.4|) . Insertion of (|5.5|) into the correlation function yields 


fsAi^o) = + do)ff'[xo) 




( 0 ), 


(5.134) 


5.4 Calculation of the critical quark mass and 
the lattice artefacts 

With the one loop expansion of /a and /p, one is now able to compute the critical 
quark mass at 1-loop order. To this end, the tree level and 1-loop coefficients of 
mi have to be computed at a given mg. From (p.l7| ) one gets the Xq dependent 
mass m{xo) at tree level by 


m 


(0) 


( 2 ^ 0 ) = 


(^0 


3o)/i"Ao) 


4/p°To) 


(5.135) 


while at 1-loop level it is given by 


m^^\xo) = 


(dl 


+ do)f^\xo) 


4^fp\xo) 


f(i) 

/p 


(a^o) 


(di 


+ do)f^\xo) 




(xo) 




d*,dofl°\xo) 

‘^fp\xo) 


(5.136) 


with and fp^ including the volume and boundary improvement terms (|5.123|) , 
(|5.127D , and ( |5.133|) . The current quark mass mi is then computed at tree and 1- 
loop level by taking m^°) and in the centre of the lattice according to (|5.18|) . 
Of course, the expansion coefficients of /a and fp depend on the bare quark mass 
at which they are computed. Thus one gets an expansion 


mi = mS°^(mo) m)"^(mo)5'o + O(5'o 


(1)^ 


(5.137) 
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where the expansion coefficients depend on the bare quark mass mo- To set up 
perturbation theory for the critical mass, we consider mo also to be a series, 

mo = m[,°^ + + 0 ( 5 ( 0 ), (5.138) 


and expand mi further as 


mi 


= m 


(0) 



+ 




I (1) 

+ mg 


d 


drui 


-m 


(0) 



9^o + 0{g^o)- 


(5.139) 


Therefore, the computation of mc^^ has to be done in two steps. First we compute 
by requiring 

(m®) = 0, (5.140) 

then we can determine mc^^ from 

mi^^ ("i®) +m^^^7-^m® = 0. (5.141) 

OTTIq 

The hrst step is easily done numerically using a discretised version of the Newton- 
Raphson method. The second step mainly amounts to expanding /a and /p up 
to 1-loop order as described in the previous subsection. Results for and mc^^ 
can be found in chapter ^ 

The lattice artefact e may now be computed at tree and 1-loop level. In 
lattice units, e is given by 


e(L) = mi(2L) — mi(L), (5.142) 

where mi is to be taken at the critical bare mass at lattice size L, which means 
that the bare mass is chosen such that mi(L) = 0. So, in order to get e(L), one 
hrst has to compute the critical mass at lattice size L and then, using this bare 
mass, compute mi on a lattice twice as large. The results for e are shown in 
chapter ^ 

The lattice artefact 

d{L) = m 2 {L) — mi{L) (5.143) 

remains to be computed. To this end, one has to expand the functions and /p 
up to 1-loop order. However, this turns out to be an easy task. It is straightfor¬ 
ward to show that and /p are the same as /a and /p with inverted background 
held. This means, in order to compute m 2 {L), one only has to compute mi and 
then repeat the computation after exchanging the boundary values C and C. 
(One has of course to take care of the 0a(a^o) in the propagators and vertices, 
which depend crucially on the background held.) The results for d can be found 
in chapter |^. 
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( 1 ) 


and fp^ 


63 


5.5 


Numerical computation of and 


Since most of the computation merely consists of matrix multiplications, the 
structure of the program is rather simple. However, due to the complicated 
formulae for the propagators and vertices, there is a big risk of errors. For 
this reason, a careful check of the results is necessary. To this end, the gauge 
parameter Aq was left arbitrary, and it was checked that the results do not depend 
on Aq. While this only gives a check on the total result, there is a possibility to 
check at least those diagrams separately that come in two types, labelled a and 
b here. Due to the symmetries of the propagators and vertices, these two types 
have to give the same results. The results using the tree level improvement 
coefficients as well as the gluon boundary term have been compared to Monte 
Carlo results obtained by Juri Rolf. Finally, two sets of programs have been 
written independently, one by Peter Weisz and one by the author of this thesis, 
and the results have been compared. 

While the general structure of the program is relatively simple, the compu¬ 
tation of the propagators is technically somewhat more involved. In contrast to 
the case of a vanishing background held, the propagators for the gluon, ghost, 
and fermion helds cannot be computed analytically. For this reason, a numerical 
calculation is required. Of course it would be possible to invert the operators JF“, 
and D by a simple inversion routine. However, on large lattices, this would 
mean inverting very large matrices, which would take a large amount of com¬ 
puter time. A more convenient way to compute the propagators is a procedure 


using Wronskian forms described in for the gluon and ghost propagators and 
in for the fermions. This method, which is briehy outlined in appendix P, was 
used in the computation of /a and /p at 1-loop order. However, even with this 
method, computer time is far from being negligible. On a 200 MHz Pentium PC, 
the computation at L = 16 took 16 hours of CPU time. Even on a 900 MHz PC, 
lattices with L > 32 could not be computed in reasonable time. Since one has 
to sum over three momentum components and two vertex times in the diagrams, 
the time needed scales asymptotically with L^. This means that larger lattices 
up to L = 64 seem to be beyond reach. 
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Chapter 6 

The renormalised quark mass 

6.1 The renormalisation constant Zp and its step 
scaling function 

6.1.1 The renormalisation constant Zp 

In chapter |^, the renormalised quark mass dehned by the PCAC relation was 
found to be related to the current quark mass via ( ^.igp . The actual renormali- 
sation is done by renormalising the axial current and density, leading to a running 
mass, since the renormalisation constants depend on the renormalisation scale. 
This is, however, not true for Za, which is only present due to the violation of 
the chiral symmetry and can be shown to be hnite and scale independent. Thus 
the renormalisation constant that actually makes the running is Zp. In order to 
compute Zp, one may use the correlation function /p. However, in /p, one does 
not only have to renormalise the pseudo scalar density, but also the boundary 
helds ( and (, such that the renormalised correlation function becomes 

(/p)R(a^o) = ZpZ'^fp{xo). (6.1) 

In order to cancel the unwelcome renormalisation factor Z|, one may use the 
correlation function /i, dehned by 

/i =-^ ^ (c'(u)75^r“C'(v)C(y)75^r“C(z)y (6.2) 

u,v,y,z ' ' 

which is renormalised by Z^ only. 


(/i)r = ^c/i- (6-3) 

This correlation function is proportional to the probability amplitude that a quark 
antiquark pair created at Xq = 0 propagates to Xq = T, a. situation which may 
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Xo= 


Xo= 



LxLxL 

Figure 6.1: The correlation function fi 


be depicted by figure pTj Like /a and /p, fi may be expanded in powers of the 
bare coupling, 

/i = f!°’ + sSft’ + 0(sS). (6,4) 

Now, Zp can be computed as 




/f (T/2) 


1 2 


h 


/<“> [MTmf 


(6.5) 


where the normalisation is chosen such that Zp is 1 at tree level. 

The renormalisation constant Zp may be expanded in powers of the bare 
coupling, 

OO 

Zp{go,L/a) = 1 + Y,Z^n\L/a)gf. ( 6 . 6 ) 

k=l 

The 1-loop coefficient Zp^ contains a logarithmic divergence. It may hence be 
decomposed into a logarithmic term, a constant term, and terms that vanish in 
the continuum limit. 


Zp^ = C-pZp{6, p) — do ln{L/a) + 0{a/L). 


(6.7) 


Here, Zp denotes the term independent of L/a. It only depends on the phase 9 
and on the ratio of the lattice extensions in time and space directions, p = T/L. 


6.1.2 The step scaling function of the running mass 

In analogy to the coupling, the running of the mass may be described by a step 
scaling function, giving the change of the running mass when changing the scale 
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L by a certain scale factor. In the following, this scale factor is chosen to be 2, 
which is a convention widely used in numerical computations. Since the running 
of the renormalised mass is completely described by the running of Zp, a suitable 
dehnition of the step scaling function at vanishing renormalised quark mass is [|75 


'£p{gl,a/L) = 


Zp{go, ‘2.LIa) 


Zp{go,L/a) 


( 6 . 8 ) 


go=9o{gR,L/a) 


with the continuum limit 


lim Sp(^^,a/L) = ap{gl). (6.9) 

a/L—>0 

As the renormalised coupling, one here has to take the Schrodinger functional 
coupling at length scale L, i.e. gp^ = g{L). Like the step scaling function a is used 

as a discretised version of the /3-function, the step scaling function up serves as 

a discretised r-function. 

The step scaling function may be expanded in perturbation theory, 

^p(5'R) a/L) = 1 + k{L/a)gp + 0{gp). (6.10) 

Since the tree level value of Zp is 1 at all scales, the 1-loop coefficient k{L/a) is 
simply given by 

k{L/a) = Z^^\2L/a) - 4^^(L/a), (6.11) 

where Z^^ihjd) and Z^^i^hjd) have to be taken at the same bare quark mass. 
In the case of a vanishing physical quark mass, this is the critical mass at length 
scale L. For the continuum limit, one gets from (|6.7| ) 

k{oc) = —(ioln(2). (6-12) 


Like in the case of the renormalised coupling, one has lattice artefacts making 
the step scaling function deviate from its continuum limit. These lattice artefacts 
can be estimated at 1-loop order by 


5k{L/a) = 


k{L/a) 

k{oo) 


1 . 


(6.13) 


In 1^, these lattice artefacts have been computed for several choices of 9 and p 


in the case of a vanishing background held. To do this calculation with a non 
vanishing background held is one of the main aims of this thesis. The size of 
the lattice artefacts is also important for the question whether it makes sense to 
compute the renormalised quark mass non-perturbatively with a non-vanishing 
background held. So far, this calculation has only be done at zero background 
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6.1.3 The 2—loop anomalous dimension 


It was already shown in subsection |2.4.3| that one may obtain the 2-loop anoma¬ 
lous dimension di in the Schrodinger functional scheme from a 1-loop computa¬ 
tion by using (|2.42|) to convert the known anomalous dimension in the MS-scheme 
to the Schrodinger functional scheme. In order to do this, one has to convert the 
coupling and the masses. For the coupling, one has 

gL), (6.14) 

where the 1-loop coefficient of Xz given by 


= 26oln(gL) - ^(cpo + Ci,iA^f). 


(6.15) 


The coefficient Ci^ has been computed in and Ci^i in |^. For 9 = ti/5 and 
p = 1, one obtains 

Ci,o = 1.25563(4), ci,i = 0.039863(2). (6.16) 

In order to get the conversion of the renormalised mass from the MS-scheme 
to the Schrodinger functional scheme, one has to use the relation between the 
renormalised mass and the subtracted bare mass. In the MS-scheme, it has been 
computed in |]^ and since then been verihed numerous times. It turns out to be 

^Ms(?) = + fi'o -^offi(ag) + 0.122282(1 )Cf -FO(5fo)|. (6.17) 

By calculating the renormalised mass in the Schrodinger functional at a non¬ 
vanishing mass, one gets |p^ 


msF(F) = mq n -hPo 


hi) 


Z^^\L/a) + 0.067886(1)Cf + O(po) 


"A 

'( 1 ) 


Inserting the results for from (b.l3|) and Zp^ from (|6.7D then yields 


msF(F) = mq n-hPo do ln(L/a) - (zp-h 0.019458(1))Cf + Oig^] 


(6.18) 


(6.19) 


Using these results for the masses in the MS scheme and in the Schrodinger 
functional scheme, one hnds the relation between both schemes. 


msF(F) = mi^(g) |l-hPo doln(gL) - (zp 4- 0.141740(2 ))Uf -4 0(5(o)|- (6.20) 

From this relation, one may easily read off the 1-loop coefficient Xm\ 

yW = do In(gL) - (zp + 0.141740(2))C'f. (6.21) 

Now, one can compute the 2-loop anomalous dimension in the Schrodinger func¬ 
tional scheme. From ( 2.42 ), one has 

df = df^ - dox^g^ + 2box^rt\ ( 6 - 22 ) 

with bo, do, and from ( p.26| ), (|2.34|) , and ( |2.35| ), and and Xm^ from (|6.15|) 
and (|6^). 
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6.2 The renormalisation constant Zj> at 1—loop 
order 


6.2.1 Preliminaries 


The expansion of /i at 1-loop order has been outlined in for the case of a 
vanishing background held. In the following, the calculation will be done for the 
non-zero background held used to dehne the coupling. 

In order to make contact to the calculation of /a and /p and thus simplify 
the computation, it is useful to rewrite /i in a slightly diherent form. To this 
end, we dehne 

^3 


K = ct 


L3 


J2{p+u{x,o)-^h{x)} 


xo=T—a 


Using this notation, fi may be written as 

/i = ^ (Ti{A-*A'}), 


(6.23) 


(6.24) 


The expansion of /i thus amounts to expanding K. This makes part of the cal¬ 
culation exceptionally easy, because all contributions in which the link variables 
only appear at tree level may be reduced to /a and /p, so these diagrams may be 
computed using the results of chapter |^. This is easily seen at tree level, where 
one has 

(6.25) 




XQ=T—a 


The clou to reducing fi to /a and /p lies in the fact that (x) does not depend 
on the spatial position x. This means that the sum over x is trivial, giving 


= F+ F(°)(x)| 


xo=T—a 


or, inserted in 


/(“) = 1 Tr 
2 




xq=T— a 

Comparing this expression to (p.30|) and (b.31|) then yields 


/!"’ = § |_/p"'(r - - /rcr - a) 

which may be decomposed into colour components. 


p(0) 


( 0 )/ 


- a) - fZ(T - a) 


1 


2 L' 


p(0) 


( 0 )/ 


(6.26) 


(6.27) 


(6.28) 


(6.29) 


The same will apply for several diagrams at 1-loop order. At tree level, /i may 
be visualised by the diagram in hgure 
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T • 


T-a 


a 




Xo=0 • 






Figure 6.2: Diagram for fi at tree level 


6.2.2 The correlation function fi at 1—loop order 

For simplicity, we return to lattice units and take a = 1 for the rest of this 
chapter. 

At 1-loop order, /i can, just like /a and /p, be written as a sum over diagrams, 
improvement terms, and the mass derivative. 


Jl — / , il,n W Cswil.V 


+5i‘’ (/ih+4/r>)+ 


Gb 


+m 


( 1 ). 


d 


dmi 




(6,30) 


The term is there because fi gets an overall factor cf from K and H. 

However, this term may be neglected, since it is cancelled by the overall factor 
cl of /p. 


For the 1-loop diagrams, one has to expand the link variables and H in ( 5.23 ) 
up to order g^. It was already stated that the diagrams containing the tree level 
values of the link variables are obviously related to the corresponding diagrams 
for /a and /p. These diagrams either contain a combination of and or of 
and H^^'>. The expression becomes particularly simple for the hrst case since 

//(o) 

is spatially constant and does not depend on the gluon fields. The gluonic 
expectation value is thus only to be taken over This expectation value 

{H^'^\x))o is again spatially constant, which means that the same argument may 
be used as in the tree level case. One thus gets for n = la, lb, 3a, 3b, 5a, 5b, 6a, 6b 


J l,n 2 


/S(r -1) - fZiT - 1) 


( 1 ) 


IP,n 


A,n\ 


(6.31) 
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T • 



6a 
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a 

C 

a 


Xo=0 

1b 

Xo=0 

2 

T 


T 


T-a 


T-a 

u 

a 

(. 

a 


Xo=0 

3b 

Xo=0 

4a 

T 


T 


T-a 


T-a 


a 


a 


Xg=0 

5a 

Xo=0 

5b 

T 


T 


T-a 

d 

T-a 


a 


a 


Xo=0 

6b 

Xo=0 

7 


Figure 6.3: Diagrams contributing to fi at 1-loop order of perturbation theory. 
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For the combination of with the situation is slightly more compli¬ 
cated. Here, the expectation value H^^\x))q does depend on x and the 

sum has a non trivial effect. When transforming the expectation value to the 
momentum representation, the sum together with the exponential factor acts as 
a delta function, leaving only the contribution with zero loop momentum. One 
thus gets for n = 2, 4a, 46, 7 


f(l) _ f (rp 

Jl,n 2 ■'P,nO 


dl) 


1) - fAi(T - 1) 


' A,n 


q=0 


(6.32) 


where [.. .]q=o means that, instead of the sum over q in (|5.102|) , (|5.105|) , (^.106|) , 
and ( b.lll|) , one has to take the term with q = 0 only. The 1-loop terms covered 
so far are represented by the diagrams in hgure |6.3i 

For the tadpole diagrams, one hnds that {H^^\x)[S^ll]p)o is independent of 
X. One thus gets 


f ( 1 ) ^ f ( 1 ) (m 

J l,n r) lp,nO 


1 ) 


1 ) 


(6.33) 


^(1) (rp 

2 J A,n\ 

for n = 8a, 8b, 9a, 9b. 

Now one is left with the diagrams containing the link variables at the upper 
boundary of the box, i.e. the link variables in (|6.23|) , at order go and Qq. The 
Qo terms can be treated in total analogy to diagrams la and lb, resulting in 
diagrams 10a and 10b of hgure They are given by 


f!!L = 


a 

OL'I 


(6.34) 


and 


HI) _ f(i) 

11,106 ~ i 1 , 10 a- 


(6.35) 


Contracting the terms of both link variables proportional to go yields diagram 11, 

di) _ A r(0) 

L3 


/tn = 71 E 4°’ E ^So(0; r - 1, T - 1)C. 


(6.36) 


Next, one has to contract the link variables at the upper boundary with Ffd) 
or, to be more precise, with the hrst order terms of the propagators and link 
variables contained in iifd) Combining them with the quark lines on the same 
side gives diagrams 12a and 12b, 

dC = 57;EE E E^H<i;“«.r-l)Tr{F+75S(0;l,r-l)75F+ 

q p SQ,tQ,uo a 

■/“S'(q; T - 1, So)V)f(q, 0, -q; Sq, h, Uo)S{0; to, 1)|, (6.37) 

fi%b = Y l>'«o)Tr|p+75F(0;l,so) 

q p sci,tci,uo a 

So,^o,Mo)<S'(q; to,T - 1)/“75P+S'(0;T - 1,1)|,(6.38) 
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13a 


13b 


Figure 6.4: Further diagrams contributing to fi at 1-loop order of perturbation 
theory. 
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while the combination with the quarks on the opposite side result in diagrams 
13a and 13b, 

M so,to,uo a 

■v;(0, 0,0; So, t.„)S(0; f„, T - l)r7oP+S(0; T - 1,1)|, (6.39) 
fl'L = E E^Ss(0;^’-l.“«)Tr{f’+7sS(0;l.r-l)75F+ 

/I so,to,uo a 

■/“E(0; T - 1, So)i;"(0, 0, 0; Sq, to, Wo)^(0; to, !)}• (6.40) 

The remaining combinations are those of the link variables at the upper boundary 
with the link variables at the lower boundary. They are depicted by the diagrams 
in hgure and can be computed via 

filL = l;^EE^So('li°.r-l)Tr{F+75S(0;l.r-l) 

q a 

■75nrs(q;r-l,l)r}, (6.41) 

= 5;^ E E "Soli; ^ 0)Tr{f>+75S(q; 1, T - 1) 

q a 

■ 7 sP+rS( 0 ;T-l,l)r}, (6.42) 

/I'.E = 

a 

■75P+/“^(0;T-l,l)r}, (6.43) 

and 

/l,15fe ^ /l,15a- (6.44) 

The only remaining diagrams are now the contractions of the hrst order terms 
of the link variables at the upper boundary with the hrst order term of the total 
action. In complete analogy to diagrams 9a and 9b one gets 

f 1,16a ~ ~2 EEE«So(0; T-l,uo)Tr{p+75^(0;l,T-l) 

fi iLQ a 

■75P+/“^(0;T- 1,1)}t;(uo), (6.45) 

and 

/l,16b ~ “/l,16a- (6.46) 

So these diagrams cancel in the sum and leave diagrams 8a and 8b as the only 
Nf dependent contributions to /i. 
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Figure 6.5: More diagrams contributing to fi at 1-loop order of perturbation 
theory. 
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Figure 6.6: Diagrams for /} y 


6.2.3 The improvement terms 

In the computation of /i, all improvement terms can be reduced to the corre¬ 
sponding terms of /a and /p. For the terms proportional to ciw and which 
are additional contributions to diagrams 6a and 6b, one may use the same ar¬ 
gument applied there, namely the independence of and {H^‘^\x))o of the 

spatial position x. One thus gets 


/i,v - 2 [/p,v(^ - 1) - fA,yiT - 1)] 


for the volume term and 


/l,Fb — - [fp,Fh(T — 1) — fA,FhiT — 1)] 


(6.47) 


(6.48) 


for the quark boundary term. These terms can be depicted by the diagrams in 
hgure |6^. 

For the boundary counter-term proportional to q , one gets several diagrams 
in total analogy to the calculation for /a and /p. For the diagrams shown in 
hgure 1^, the same argument applies as above, yielding 


/l,Gb — - [/p,Gb(^ “ 1) “ fA,GhiT — 1)] . 


(6.49) 


Apart from these diagrams, one gets the contributions shown in hgure |6]^. How¬ 
ever, like for /a and /p, these contributions cancel pairwise in the sum and may 
thus be ignored. So the improvement terms for fi can easily be computed us¬ 
ing (^T^ , (^TTH) , and ( Pr^) . 

Also for /i, we will need the derivative of the tree level coefficient with respect 
to the bare mass. Diherentiating (|6.28|) gives 


d 


/!”> = 


1 d 


h"’(r-i)-/r(r-i) 


( 0 ), 


dmo 2 dmo l 

do) _ 1 .(0) 


(6.50) 


where /p and can be taken from ( 5.92|) . 
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Im3a 


Im3b 


Figure 6.7: Diagrams contributing to /{ Qb 


6.2.4 Computation of Zp at 1—loop order 

The renormalisation constant Zp at 1-loop order is obtained by expanding fi 
and fp in (|675|) . One thus gets 


(1) _ /d 4‘’(r/2) 


Zp — 


2^(0) f^->^T/2) 


( 0 ), 


(6.51) 


In order to get Zp at vanishing renormalised quark mass, fi and fp, which depend 
on the bare quark mass, have to be computed at the critical quark mass rric- In 
perturbation theory, one has to use the expansion 


nrir = 


+ m^^^gl + 0{g\ 


(6.52) 


where the computation of m'c'^ and m'c'^ has been outlined in chapter ^ Taylor 
expansion then yields 

= /i(m®) + + o(^^), (6.53) 

which leads to the 1-loop coefficient 


(6.54) 
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Figure 6.8: More diagrams for Qt,. These contributions cancel in the sum. 
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and analogously for /p. 

Furthermore, and are meant to include all improvement terms. Then 
the results obtained in the calculation should converge to their continuum limits 
at a rate proportional to a?. 

By calculating all the diagrams, one is now able to compute the renormali¬ 
sation constant Zp at different lattice sizes. In the case of 0 = vr/S and p = 1, 
the coefficient defined in (|6.7| ) could be obtained by extrapolation to the con¬ 
tinuum limit, using the extrapolation procedure described in appendix p. From 
this calculation, one automatically gets the 2-loop anomalous dimension in the 
Schrodinger functional scheme by (|6.22|) . This means that the 2-loop anomalous 
dimension could be obtained by a 1-loop calculation, since the 2-loop anomalous 
dimension in the MS-scheme is already known. 

Setting mi(L) = 0 and computing Zp{L) and Zp{2L) at this mass gives the 
1-loop coefficient k{L) of the step scaling function. One is thus able to compute 
the size of the discretisation error (|6.13| ). This computation was done in the case 
oi 6 = n/5 and p = 1. For some smaller lattices, we did the calculation also at 
6 = 0 and 6 = 0.5 and at p = 2, in order to get a better comparison with the 
results obtained in [^, which were computed at these values of 9 and p. 


Since has not been computed with a non-vanishing background field 
before, there are no results to compare the results obtained in this calculation 
to. For this reason, careful checks of the results are necessary. For /i, all checks 
described for /a and fp are applicable. Moreover, it was checked that the results 
both for the sum of the diagrams and for the improvement terms are invariant 
when interchanging the boundary conditions in time direction, as they should be. 

All results can be found in chapter |^. 
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Chapter 7 
Results 


The following results have been computed at vanishing renormalised mass, which 
means setting mi = 0 or mo = rric- The results depend, of course, on the phase 
angle 9 and the ratio p = T/L. Unless stated otherwise, all results presented here 
are for 9 = ir/b and p = 1. 


7.1 The current quark mass 


7.1.1 The critical quark mass 

At tree level, the critical quark mass is independent of the flavour number, while 
at 1-loop level, there is a linear dependence on iVf due to the closed quark loops 
in diagrams 8a and 8b in figure p.4|. The 1-loop coefficient may hence be written 


as 


■nrS}^ = m^o^ + m^i^iVf. 


(7.1) 


The results for the tree level coefficient m, 

di' 


and the 1-loop coefficients rri'^Q and 
can be found in table O- The critical mass at 1-loop level is also shown 
in figure [f.l| for different ffavour numbers. The unusual case Nf = —2 refers to 
so called bermions, which were invented because Monte Carlo simulations are 
considerably cheaper in this model than in full QCD. The original idea was to 
extrapolate from negative values to iVf = 2 |^, , which is, however, problem¬ 


atical. Nowadays, bermions are more used as a toy model |^, 

A clear dependence of the 1-loop critical quark mass on the flavour number 
can only be seen on very small lattices. For larger L/a, the 1-loop coefficient can 


be seen to be practically constant at its continuum limit value, which is |^, 


am(^) = -0.2700753495(2). 


(7.2) 
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Results 


L/a 

(0) 

amc 

(1) 

aKo 

(1) 

am^i 

4 

-0.0015131 

-0.2653473 

0.0024231 

5 

-0.0016969 

-0.2675073 

0.0006225 

6 

-0.0006384 

-0.2699739 

0.0002269 

7 

-0.0005761 

-0.2700126 

0.0000978 

8 

-0.0003209 

-0.2700683 

0.0000522 

9 

-0.0002753 

-0.2700656 

0.0000300 

10 

-0.0001835 

-0.2700733 

0.0000192 

11 

-0.0001561 

-0.2700708 

0.0000129 

12 

-0.0001145 

-0.2700730 

0.0000091 

13 

-0.0000979 

-0.2700720 

0.0000066 

14 

-0.0000761 

-0.2700731 

0.0000050 

15 

-0.0000657 

-0.2700729 

0.0000039 

16 

-0.0000531 

-0.2700736 

0.0000030 

17 

-0.0000463 

-0.2700736 

0.0000024 

18 

-0.0000385 

-0.2700741 

0.0000020 

19 

-0.0000339 

-0.2700742 

0.0000016 

20 

-0.0000288 

-0.2700745 

0.0000013 

21 

-0.0000256 

-0.2700746 

0.0000011 

22 

-0.0000221 

-0.2700749 

0.0000009 

23 

-0.0000198 

-0.2700750 

0.0000008 

24 

-0.0000173 

-0.2700751 

0.0000007 

25 

-0.0000156 

-0.2700752 

0.0000006 

26 

-0.0000138 

-0.2700753 

0.0000005 

27 

-0.0000125 

-0.2700754 

0.0000005 

28 

-0.0000112 

-0.2700755 

0.0000004 

29 

-0.0000102 

-0.2700755 

0.0000004 

30 

-0.0000092 

-0.2700756 

0.0000003 

31 

-0.0000084 

-0.2700756 

0.0000003 

32 

-0.0000076 

-0.2700756 

0.0000003 


Table 7.1: Perturbative results for the critical quark mass me up to 1-loop order 
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Figure 7.1: The critical quark mass me at 1-loop order 
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Results 


7.1.2 Lattice artefacts of the current mass up to 1—loop 
order 


Using the expansion of /a and /p, it is possible to compute the lattice artefacts 
d and e dehned in (|5.23|) and (|5.24|) . Here again, the 1-loop coefficients get a 
linear TVf dependence from diagrams 8a and 8b in hgure |5.4| and may hence be 
written as 


= d, 


( 1 ) 


0 +d?Nf, 

eW = 


(7.3) 

(7.4) 


The results at tree and 1-loop level can be found in table [7^ , the 1-loop coefh- 
cients are also shown in hgure |7.2| for different values of Nf. They are obviously 
small and converge quickly to zero. The expected behaviour in the improved 
theory is a convergence at rate proportional to (a/L)^. Figure |^, where the 
1-loop coefficients are plotted with respect to (a/L)^, does however show that, 
for iVf 7 ^ 0, they are dominated by terms of higher order. In |^|, the perturbative 
results up to 1-loop order are compared to results obtained by Monte Carlo simu¬ 
lations. While the values for e are approximated by perturbation theory remark¬ 
ably well, the perturbative values for d are signihcantly larger than the simulation 
results, especially at large values of the renormalised coupling {g^{L) = 3.334). 

An important question in this context is the dependence of the results on the 
improvement term proportional to which was missing in previous publica¬ 
tions. As already stated in chapter |^, this term is absent in zero background held 
calculations and irrelevant for the extrapolation both of Csw and It has, 

however, to be taken into account in the results for d{L/a) and e{L/a) published 
^|. There, the perturbative results including tree and 1-loop level are plot- 


m 


ted for different values of the renormalised coupling g^. The results for the same 
values of g^ are shown in hgure [f.3|, where the symbols denote the results with 
full improvement, while the lines show the results without the c^^^-term, i.e. the 
values presented in ||^. The picture shows that the additional term only gives 
considerable changes in the case g^ = 3.334 for W = 0 and iVf = —2. Since, for 
this coupling, only the results for N{ = 2 are plotted in 
the plots presented there are almost invisible. 


271, the corrections to 


7.2 Lattice artefacts of the step scaling function 
of the coupling 


One of the reasons to compute the critical quark mass up to 1-loop order was 
to get the correct discretisation errors for the coupling up to 2-loop order. To 
this end, one has to use ([4.45| )-( p:.4f]|) with the current mass properly set to zero, 
mi{L) = 0. Since the critical mass converges quickly to its continuum limit. 
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L/a 


Wi" 

LdP 


LeH’ 

LeP 

4 

0.1123 

0.00752 

0.01929 

-0.004981 

0.01931 

0.00990 

5 

0.0522 

0.00535 

0.00653 

-0.007783 

0.01252 

0.00309 

6 

0.0258 

0.00378 

0.00348 

-0.003204 

0.00032 

0.00133 

7 

0.0140 

0.00320 

0.00172 

-0.003550 

0.00011 

0.00066 

8 

0.0078 

0.00277 

0.00112 

-0.002165 

-0.00015 

0.00040 

9 

0.0044 

0.00224 

0.00071 

-0.002149 

-0.00012 

0.00025 

10 

0.0025 

0.00180 

0.00053 

-0.001558 

-0.00013 

0.00018 

11 

0.0014 

0.00149 

0.00038 

-0.001483 

-0.00009 

0.00013 

12 

0.0007 

0.00122 

0.00030 

-0.001172 

-0.00008 

0.00010 

13 

0.0002 

0.00104 

0.00023 

-0.001098 

-0.00006 

0.00008 

14 

-0.0001 

0.00087 

0.00019 

-0.000913 

-0.00005 

0.00006 

15 

-0.0002 

0.00076 

0.00016 

-0.000850 

-0.00004 

0.00005 

16 

-0.0003 

0.00065 

0.00014 

-0.000730 

-0.00003 

0.00004 

17 

-0.0004 

0.00057 

0.00011 




18 

-0.0005 

0.00050 

0.00010 




19 

-0.0005 

0.00045 

0.00008 




20 

-0.0005 

0.00040 

0.00007 




21 

-0.0005 

0.00036 

0.00007 




22 

-0.0005 

0.00032 

0.00006 




23 

-0.0005 

0.00029 

0.00005 




24 

-0.0005 

0.00026 

0.00005 




25 

-0.0005 

0.00024 

0.00004 




26 

-0.0004 

0.00022 

0.00004 




27 

-0.0004 

0.00020 

0.00003 




28 

-0.0004 

0.00019 

0.00003 




29 

-0.0004 

0.00017 

0.00003 




30 

-0.0004 

0.00016 

0.00003 




31 

-0.0004 

0.00015 

0.00002 




32 

-0.0003 

0.00014 

0.00002 





Table 7.2: The lattice artefacts d = m 2 {L) — mi{L) and e = mi(2L) — mi{L) up 
to 1-loop order 














86 


Results 



Figure 7.2: The lattice artefacts d = m 2 {L/a) — mi{L/a) and e = mi{2L/a) — 
mi{L/a) at 1-loop order 
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(a/L)" (a/L)2 


Figure 7.3: The lattice artefacts e = mi{2L/a) — mi{L/a) and d = m 2 {L/a) — 
mi{L/a) up to 1-loop order at different values of the renormalised coupling g"^. 
The lines show the values without the c[^^-term, while in the symbols, it is in¬ 
cluded. 
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the discretisation errors 6 ij, shown in table |7.3| are expected to differ from the 
estimation in |Q, which was obtained using the continuum limit value of the 
critical mass, only for small lattices. Since <520 does not depend on the critical 
quark mass, the 2-loop coefficient will remain unaltered for iVf = 0 when using the 
correct critical mass instead of its continuum limit. From the size of the critical 
mass (figure |^), one expects a small deviation for iVf = —2atL/a = 6, while 
for Nf = 2, larger deviations are expected at L/a = 4 and L/a = 6 . Figure M 
shows that this indeed is the case. The lines are the estimates from [^, while 
the symbols denote the values obtained with the proper mc{L/a). In any case, 
the 2-loop coefficient 62 is small and does vanish at a rate proportional to {a/L)‘^, 
as one should expect in the improved theory. Especially for iVf = 2, it does get 
even smaller by using the correct critical mass. 


L/a 

^10 

5ii 

<^20 

<^21 

S 22 

4 

-0.01033 

0.00002 

-0.001588 

-0.000630 

0.000712 

5 

-0.00625 

-0.00014 

-0.000872 

-0.000460 

0.000407 

6 

-0.00394 

-0.00014 

-0.000550 

-0.000341 

0.000200 

7 

-0.00268 

-0.00014 

-0.000376 

-0.000211 

0.000103 

8 

-0.00194 

-0.00011 

-0.000268 

-0.000138 

0.000059 

9 

-0.00148 

-0.00009 

-0.000196 

-0.000098 

0.000038 

10 

-0.00117 

-0.00007 

-0.000147 

-0.000073 

0.000027 

11 

-0.00095 

-0.00006 

-0.000111 

-0.000058 

0.000020 

12 

-0.00079 

-0.00005 

-0.000085 

-0.000047 

0.000016 


Table 7.3: The discretisation errors Sij up to 2-loop order 


7.3 The renormalised quark mass 

7.3.1 The renormalisation constant Zp up to 1—loop order 

In order to obtain the finite part Zp of Zp defined in ( |6.7|) , one has to extrapolate 
the quantity 

Yp{L/a) = ^ (4'^(L/a) + doln(T/a)) (7.5) 

to the continuum limit. Due to the iVf dependent tadpole graphs, one has a linear 
Nf dependence of Yp, 

Yp = Ypo + NfYpf, (7.6) 

and thus the same dependence for Zp, 


^po T -^f'^pi- 


(7,7) 
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Figure 7.4: The 2-loop discretisation error 62 - The lines show the values for 62 if 
the continuum limit of the critical mass is used. 
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L/a 

hpo 

hpi 

4 

-0.13990 

0.025124 

6 

-0.14008 

0.009000 

8 

-0.13757 

0.005284 

10 

-0.13590 

0.004070 

12 

-0.13476 

0.003539 

14 

-0.13395 

0.003254 

16 

-0.13336 

0.003081 

18 

-0.13291 

0.002966 

20 

-0.13256 

0.002886 

22 

-0.13228 

0.002828 

24 

-0.13205 

0.002785 

26 

-0.13187 

0.002751 

28 

-0.13172 

0.002724 

30 

-0.13159 

0.002703 

32 

-0.13148 

0.002686 


Table 7.4: The Nf independent and Nf dependent part ofYp(L/a) 


The results for different lattice sizes are shown in table [7.4| . The behaviour of Yp 
is shown in hgure for different values of Np 

For the extrapolation, the procedure described in appendix y may be used. 
Following [^, one expects an asymptotic expansion 


lp(L/a) = ao + K + h ln(a/L)] (a/L)’^ (7.8) 

k=l 


to hold, where the terms linear in {a/L) are cancelled by improvement. Taking 
the functions 1, {a/L)"^, and (a/L)^ ln(a/L) into the £t and using {a/L)^ and 
{a/L)^ ln{a/L) for the estimation of the systematic error yields the results 


Zpo = -0.13044(2), Zpi = 0.002565(10). (7.9) 

Since Yp depends on the improvement coefficients ciw\ c^\ and via the critical 
mass also on the errors of these coefficients have to be taken into account. 
However, varying the improvement coefficients within their error ranges shows 
that the change in Zpo and Zpi is smaller than the errors quoted in ([7.9|). The 
error due to the uncertainty in the improvement coefficients is thus negligible. 

Taking (a/L)^ and (a/L)^ ln(a/L) into the £t and using the functions (a/L)'^ 
and (a/L)^ ln(a/L) for the error estimation does not change the picture signifi¬ 
cantly. While one gets smaller systematic errors in the extrapolation, the errors 
of the improvement coefficients become more important. Taking everything into 
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Figure 7.5: The finite part Yp of the renormalisation constant Zp at 1-loop order 
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account, the errors quoted in (|7.9|) seem to be reasonably conservative. Including 
{a/L) and (a/L)ln(a/L) in the fit shows that the coefficients of these functions 
are compatible with zero, as one should expect in the 0(a) improved theory. 

Inserting Zp into (|6.21| ) and using ( |6.22|) , one now gets the 2-loop anomalous 
dimension . Due to the linear dependence of Zp on iVf, the anomalous dimen¬ 
sion will get an N{ dependence up to order N^, in contrast to order N{ in the 
case of a vanishing background field. Now, becomes 


df/do = 0.16508(8) - 0.00724(4)iVf + 0.000570(2)iV^. 


(7.10) 


This result for 6 = tt/5 and p = 1 is small, but larger than the anomalous 
dimension with a vanishing background field computed in |^, which is df^/do = 
0.1251 + O.OOAGNi for 0 = 0 and df/do = 0.0271 + 0.0105iVf for 6 = 0.5. 


7.3.2 Lattice artefacts of Ep up to 1—loop order 

By computing the renormalisation constant Zp both at L/a and 2L/a with 
mi{L/a) = 0, one obtains the one loop coefficient k{L/a) of the step scaling 
function Sp, 

k{L/a) = Z^p\2L/d) - zf{L/d), (7.11) 

and its deviation from the continuum limit, 

4(i/«) = - 1, (7.12) 

with fc(oo) = —do ln(2). In contrast to the zero background field case, the tadpole 
graphs give a dependence on Nf, 


5k — 5ko + 5ki-5lf, 


(7,13) 


where the coefficients 5ko and 5^1 (not to be confused with the 6ij from section [f.2|) 
are expected to vanish at a rate proportional to a?. The results for 6ko and Ski 
can be found in table |^. In hgure |7.6| , Sk is shown for several flavour numbers. 
Compared to the results obtained with a vanishing background field in [^, these 
discretisation errors are very large, especially for L/a = 4 and L/a = 6. However, 
for iVf = 2 and L/a> 8, the error seems to be reasonably small. 

In order to get a direct impression of the change of the discretisation errors 
due to the presence of the background field, one should compare the values with 
and without background field at the same values of 6 and p. This is done in 
table and in figure where the black lines represent the iVf independent 
results with a vanishing background field from |^, while the symbols denote 


the results obtained with the non-vanishing background field for several flavour 
numbers. For iVf = 0, the discretisation errors with a non-vanishing background 
field are smaller than those with zero background field at 6 = 0, while at 6 = 0.5, 
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L/a 

6ko 

6ki 

4 

-0.3258 

0.6248 

6 

-0.1998 

0.1916 

8 

-0.1548 

0.0790 

10 

-0.1229 

0.0428 

12 

-0.0998 

0.0274 

14 

-0.0826 

0.0193 

16 

-0.0696 

0.0144 


Table 7.5: The Nf independent and N{ dependent part of 6k{L/a) for 6 = 7r/5 
and p = 1 



Figure 7.6: The 1-loop order discretisation error 6k of the step scaling function 
Sp for 9 = tt/ 5 and p = I 
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Figure 7.7: The 1-loop order discretisation error 6k of the step scaling function 
Sp for various values of 9 and p at Nf = 2 (o), Nf = 0 (\/), and Nf = —2 (u). 
The solid line (—) shows the values for a vanishing background held from |3^/. 
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L/a 

(5fco e=o 

0=0 

0=O,V=1 

0 = 0.5 

0 = 0,5 

0 = 0 . 5 ,V =1 




P = 

1 



4 

0.1425 

0.2688 

0.2040 

-0.1697 

0.5765 

0.2136 

6 

0.0538 

0.1046 

-0.0121 

-0.1167 

0.1891 

0.0208 

8 

0.0164 

0.0442 

-0.0253 

-0.0983 

0.0791 

-0.0026 

10 

0.0033 

0.0231 

-0.0215 

-0.0815 

0.0425 

-0.0062 

12 

-0.0019 

0.0144 

-0.0171 

-0.0680 

0.0270 

-0.0064 

14 

-0.0040 

0.0100 

-0.0137 

-0.0573 

0.0190 

-0.0058 

16 

-0.0049 

0.0074 

-0.0111 

-0.0489 

0.0141 

-0.0052 




P = 

2 



4 

0.0963 

0.5651 

-0.3084 

-0.5116 

3.1960 

-0.2456 

6 

-0.0020 

0.1768 

-0.2292 

-0.3020 

1.1311 

-0.1499 

8 

-0.0144 

0.0880 

-0.1449 

-0.2007 

0.4744 

-0.0893 

10 

-0.0155 

0.0546 

-0.0974 

-0.1472 

0.2502 

-0.0584 

12 

-0.0142 

0.0377 

-0.0696 

-0.1137 

0.1567 

-0.0412 


Table 7.6: The iVf independent and Nf dependent part of Sk{L/a) for various 


values of 6 and p compared to the zero background held {V = 1) values from |p5 


they are larger. For iVf ^ 0, the discretisation errors become drastically larger due 
to the large iVf-term, which is not present in the case of a vanishing background 
held. This ehect is particularly large at 6^ = 0.5 and p = 2. 


7.4 Numerical checks of the results 


In order to check the programs, they were used to compute the tree level and 
1-loop coefhcients at mo = 0 and L/a = 6, where the results can be compared to 
non perturbative results obtained in Monte Carlo simulations by Juri Rolf. The 
simulations were done at small bare couplings ranging from 0.075 to 0.6 using 
the tree level values for the improvement coefhcients. The results are shown in 
hgure 1^. The tree level values of /p(T/2) and fi are obtained as the constant 
parts in polynomial hts up to order g^. (The results do, of course, depend on the 
ht ansatz. Using a polynomial of the form oq + aig^ + 02(70 turned out to give 
results with reasonable errors.) Fitting the polynomials to the hlled symbols in 
hgure yields 

f^\T/2) = 0.8586(11), = 0.1312(3), (7.14) 


compared to the perturbative results 

f'i\T/2) = 0.8587, 


= 0.1310. 


(7.15) 
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Figure 7.8: Monte Carlo results for /p(T/2) and fi at L/a = 6. The error bars are 
smaller than the symbols. The lines are polynomial fits up to order Qq including 
the hlled symbols, while the open symbols were disregarded. 


Obviously, these results agree perfectly well. 

For the one loop coefficients, one could in principle use the coefficient of 
obtained in the £t. However, assuming the tree level values of fp and /i to be 
correct, one achieves more precise results doing the £t on {fp{T/2) — fp\T/2))/gQ 
and (/i — /i°^)/5'o) where the 1-loop coefficients get the constant parts in the fit. 
These functions are shown in hgure [7.9| both for tree level improvement and for 

= 2, which was used to check the c^-term. In this case, the results for bigger 
gQ seem to be dominated by terms of higher order than g^. 

For c[^'^ = 0, the extrapolation yields 

/^i)(T/2) =-0.765(11), =-0.188(3), (7.16) 


compared to 


/W(T/2) = -0.760, 


/P = -0.189. 


(7,17) 
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Figure 7.9: Monte Carlo results for {fp{T/2) — fp\T/2)) / and (/i —/i°^)/fi'o 

L/a = 6. The error bars are smaller than the symbols. The lines are polynomial 
hts up to order including the filled symbols. 

obtained by perturbation theory, while for = 2 one gets 

/W(T/2) =-0.480(10), =-0.186(3), (7.18) 

compared to the perturbative results 

fi^\T/2) = -0.466, = -0.183. (7.19) 

These results seem to agree well enough to state that the perturbative results 
could be conhrmed by the simulations. 
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Chapter 8 
Summary 


In this thesis, a part of the programme of the ALPHA collaboration to renormalise 
QCD has been presented. 

The basic ideas of the Schrodinger functional as a method to compute renor- 
malised quantities have been outlined and it has been explained how this scheme 
avoids computations on very large lattices by using the length of the space time 
box as the renormalisation scale. Furthermore, Symanzik’s improvement pro¬ 
gramme has been introduced as a method to reduce discretisation errors to O(a^). 

The main purpose of the Schrodinger functional method, the renormalisation 
of parameters, has been studied in detail. The definition of the renormalised 
coupling has been outlined, and it was shown how a renormalised quark mass 
can be introduced via the PCAC relation. 

In order to do several 1-loop calculations, perturbation theory in the Schro¬ 
dinger functional has been discussed. Perturbation theory has then been used to 
do a variety of calculations with a non-vanishing background field, which so far 
had only been done in the zero background field case. In particular, the following 
calculations have been done: 

• The correlation functions /a, /p, and /i have been expanded up to 1-loop 
order. 

• In these expansions, it was found that, with a non-vanishing background 
field, an improvement term proportional to arises, which was missing 
in previous publications. 

• The critical quark mass has been computed up to 1-loop order and was used 
to calculate the discretisation errors of the coupling’s step scaling function 
up to 2-loop order. 

• Several discretisation errors of the current quark mass have been calculated. 

• The renormalisation constant Zp has been computed at 1-loop order, giving 
the 2-loop anomalous dimension. 
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• The discretisation error of the step scaling function of the quark mass has 
been computed at 1-loop order. 

These computations turned out to be time consuming, so that lattices larger 
than L/a = 32 could not be reached. Up to L/a = 32, calculations could however 
be done in reasonable time. In these calculations, the following results have been 
obtained: 

• The critical quark mass was shown to converge quickly towards its contin¬ 
uum limit. Considerable deviations from this limit were only observed for 
very small lattices. 

• The discretisation errors of the coupling’s step scaling function turned out 
to be small and close to those computed with the continuum limit of the 
critical quark mass. Where the results using the correct critical mass deviate 
from the previous results, they were shown to make the discretisation errors 
smaller. 

• The discretisation errors of the current quark mass turned out to be small. 

• The effect of the previously missing 4^^-term on results published so far 
was shown to be negligible. 

• The 2-loop anomalous dimension was shown to be small but larger than in 
the case of a vanishing background held. 

• The discretisation error of the step scaling function of the renormalised 
quark mass turned out to be much larger than with a vanishing background 
held for iVf ^ 0, while for Nf = 0, it is smaller or larger, depending on 6. 
For the case of 0 = vr/S, p = 1, and N{ = 2, it turned however out to be 
reasonably small for L/a > 8. 

• The results obtained perturbatively could be shown to agree well with non- 
perturbative results from Monte Carlo simulations in the small coupling 
limit on small lattices (L/a = 6). 

With these calculations, the perturbative treatment of the quark mass renor- 
malisation up to 1-loop order at vanishing renormalised mass is completed. The 
programs are, however, written such that they can be used for calculations at 
diherent quark masses without too much ehort. 

What has not yet been done are non-perturbative calculations of the renor- 
malised mass with a non-vanishing background held. However, due to the large 
discretisation errors found in this thesis for small lattices, safe extrapolations to 
the continuum limit could turn out to be difhcult. 



Appendix A 

Notations and conventions 


A.l The Dirac matrices 


The Dirac matrices in Euclidean space may be obtained from those in Minkowski 
space by 


^Euclidean 
^ Euclidean 

7i,2,3 


^ Minkowski 

7o 

Minkowski 

-^7i,2,3 


giving 


7o = 


71,2,3 — 


with the Pauli matrices 


CTi = 


0 1 
1 0 


era = 


1 0 
0 -1 y ’ 

0 -*cri,2,3 

* 0 - 1 , 2,3 0 


0 -i 


1 0 

I 0 y ’ “ I 0 -1 


In this representation, 75 = 71727374 is given by 

0 -1 


75 = 


-1 0 


The projectors P+ = ^(1 + 70 ) and P_ = |(1 — 70 ) become 


P^ = 


1 0 
0 0 


p_ = 


0 0 
0 1 


(A.l) 

(A. 2 ) 

(A.3) 

(A.4) 

(A.5) 

(A. 6 ) 

(A.7) 
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A.2 The basis of the Lie algebra su(3 


One possible choice for a basis of su{3) are the well known Gell-Mann matri¬ 
ces For the calculations in this thesis, it is however more convenient to use 


the same basis as in [^. To this end, we introduce colour matrices Xa, which 
coincide with the Gell-Mann matrices Xa in all cases except for the two diagonal 
matrices, for which we choose 


^3 — “ 2 "^^ ^ 
^8 = ^“-^3 + 


Then, the matrices are 

Ai = 

A 4 = 

\ = 

^3 = 

After normalisation. 



Ao — 


Ak — 


A? — 



r 

0 

0 ^ 

\ -11 

r 2 

0 



1 

0 

I- 

0 

-1 


\ 0 

0 

-1 } 

^ 0 

0 

-1 / 


Ta — wtAq, 

2i 

these matrices may be used to dehne a new basis which is given by 


/' = A = ^(Ti + *T2), 

/- = (7+ = + iTs), 


7" = r_ = ^(Ti - ffj), 
P = U. = -L[T, - iT,), 


(A.8) 

(A.9) 


(A.IO) 


(All) 


7® = n = + ^T^), 


7' = r_ = ^(Te - iTi), (A12) 


for the non-diagonal matrices and 
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a 

Ca 


1 

1 (cos 27 + cosy) 

—il(sin 2 y + siny) 

3 

cosy 

0 

4 

|(cos 2 y + cosy) 

—z|(sin 2 y + siny) 

6 

cosy 

0 

8 

i( 2 cos 2 y + cosy) 

0 


Table A.l: Ca and Sa for the gauge group SU(3). The other coefficients are 
C 2 = Cl, C, = C,, Cr = Cq, S 2 = -Si, S, = -S^, and S 7 = -S^. 


for the diagonal ones. For this basis, one has 

= -C, (A.14) 

where 1 = 2, 4 = 5, 6 = 7, and vice versa. For the diagonal matrices, one has 
3 = 3 and 8 = 8 . The normalisation is chosen such that 

Tr (/“/') =-^5b,. (A.15) 


A.3 The background field 


In the chosen basis, the constant colour electric background field S is proportional 
to J®. It can be written as 


with 


/ 2 0 0 \ 

T = —( 0 —1 0 j = —'yVsXs, 
Vo 0 -ij 


(A.16) 


(A.17) 


A special feature of the chosen basis is that the star operation defined by 


M*G 


{MG + GM^) 


1 

N 


tr(MG + GM^) 


(A.18) 


acts diagonally. 


coshGofc^/“ = CaG, (A. 19) 

sinhGofc^r = SaG. (A.20) 


The values of Ca and Sa can be found in table |A-1 
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a 

Mxo) 


Ra 

1 

-37x0 + 

- -) 

S'’ 

COS 1 

3 

0 


cos 7 

4 

-37x0 + 7 (d[i + ^] 

2n \ 

3 ) 

cos 1 

6 

l(2n, - 1) 


cos 7 

8 

0 


i(2 cos 27 -|- cos 7) 


Table A.2: 0a(2^o) and Ra for the gauge group SU{3). The other coefficients are 
02 = ~01) 05 = ~04) 07 = “06) -^2 = -^1; -^5 = -^4) and -Ry = i?6- 


The covariant derivative may be decomposed in the basis J“, giving 


(Dt/)W = 5] k(io)r(i^+ *)-/“ win, 


(Dn)w = 5;[r(i)-ww)7“(2^-h 

a 

where 12a(a^o) is a phase factor and can thus be written as 




(A.21) 

(A.22) 

(A.23) 


The values of 0a are shown in table Another quantity useful for the calcu¬ 
lation of the gluon propagator is 


Ra = {Ca - 


(A.24) 


Also the values of Ra can be found in table A..2. 











Appendix B 

The numerical construction of 
the propagators 


B.l The gluon and ghost propagators 


In order to obtain the gluon and ghost propagators, one has to invert the operators 
and from (|5.51|) , (|5.52|) , (|5.53|) , ( b-54|) , and (|5.58|) . A very convenient way 
to do this is described in |^| and will be briefly outlined. 

For the gluon propagator, one has to solve the difference equation 


EE 

yo ^ 


This is particularly simple for a = 3,8 with p = 0. In this case, one simply has 




A(1 + min(xo, 2 /o)) 
;A(min(xo,|/o) - 

o“ 


if /i = = 0, 

ii H = p = k, 
if yU 7 ^ z/. 


(B.2) 


In all other cases, the difference equation has to be solved numerically. To simplify 
the notation, we denote the propagator at hxed values for p and a by 


Dl(p; xo, m) = Q„Axo, m)- 


(B.3) 


The operator may be decomposed according to 

^ 0 ) l/o) T (^o)*5x0, 1/0 T 1)^. 

Using this notation, the difference equation (|B.1|) becomes 


(B.4) 


*^/icr(^o)Qcri ^(^0 I? Z/o) ^Vo) 

(^0 l)Q(j//(^0 I 5 2 /o) ^fii'SxQ^yQ' 


(B.5) 
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The numerical construction of the propagators 


In order to solve this equation, we now construct two solutions and 

of the homogeneous equation by two step recursion forward and backward in time. 
The forward solution 'ijjf starts from 


V'o.(-l) = V'o.(O) 


(B.6) 

V'L(o) 

= 0, 

(B.7) 

V'L(i) 

^ku) 

(B.8) 

while the backward solution 'ip^ starts from 



4 .{T) = -1) 

= 

(B.9) 

^liT) 

= 0, 

(B.IO) 


^ku- 

(B.ll) 


(|B.5|) expresses the fact that Q is the right-inverse of expressed in terms of 
A and B. Then, Q is also the left-inverse, leading to the difference equation 


Qyo + l).Ayo-(|/o) Q/x(t(XO; yo)Bau{yo') 

“fQ/itr(^0) l/O 1) ^iJ,u^xo,yo- 

The homogeneous solutions of ( [B.12| ) are also given by and Assuming the 
validity of ( p.5|) and (p.l2|) for xq ^ y^ and demanding the symmetry 

Q yuixQ^yo) Quyiyo^ Xq) (B.13) 


leads to 


Qyu{,XQ, ?/o) 


for Xq < T/q 

^Li.Xo)y^~lpjix{yo) for Xq > 2/0 


(B.14) 


We now demand that the matrix W makes this dehnition consistent and also for 
2^0 = 2/0 solves (|B.5|) and (|B.12|) . These requirements can be shown to imply 


= '04(Xo)Xa(2:o)V'L(^O + 1) - ¥aa{XQ + 1)Aa(^^o)V'L(^o)• (B.15) 


This is a Wronskian form which is easy to compute. Using ( [B.5| ) and the fact 
that 'ipl^{xo)Aaxixo)'ip\j^{xQ + 1) is symmetric in y and W can easily be shown 
to be independent of Xq. A particularly simple choice is Xq = 0, which gives after 
some algebra 


Wo. = ^o<.(0)C(l) + (Boo(0)+Ao(-l)X(0), (B.16) 

Wfc. = -A,k{mU^)- (B-17) 


Now one only has to insert W into (|B.14|) to get the propagator. 

The computation of the ghost propagator proceeds along the same lines. Due 
to the absence of the spin matrix structure, it is, however, much simpler. 
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B.2 The quark propagator 


The quark propagator is computed by carrying over the methods used for the 
gluon and ghost fields. The propagator S is the inverse of the Dirac-Wilson 
operator D from ( b-61[) , which means one has to solve the equation 

Eo(p ]Xo,yo)S{^]yo,Zo) = (B.18) 

yo 

with the boundary conditions 

P+^(p;xo,2/o)Lg=o = ^->^(P;^o,l/o)Lo=r = (^-19) 

S{p;xo,yo)P-\y^^o = S{p;xo,yo)P+\y^^T = ^^ (B-20) 

which are obtained by setting the boundary quark fields to zero. 

In analogy to the gluon case, one first constructs solutions of the homogeneous 
equation by forward and backward recursion. This means one solves 


^(P’ ^0’ yoWip, yo) = 0 

(B. 21 ) 

yo 

^ P>(p; Xo, 2/o)V'^(P, 2 / 0 ) = 0, 

(B.22) 

yo 

0 

0 

II 

0 

H 

0 

(B.23) 

P.'iP\p,Xo)l^^^ = 0. 

(B.24) 


and 

starting from 
and 


Here, only one starting value per solution is required, since one only has to solve 
a first order difference equation. This is easily seen by defining 

F^/\xo) = xo) + P+V'^/'(P, Xo-1). (B.25) 

Writing down the Dirac-Wilson operator in the form of (|5.61|) , 

D{p-,xo,yo) = +P(p+,xo)5xo,yo -P+^xo-i,yo, (B.26) 

then gives 

[P(p+, xo)P+ - P-] F^/\xo + 1) + [P(p+, xo)P. - P+] F^^\xo) = 0, (B.27) 

which is solvable by one-step recursion. 

Having computed the solutions and one can construct the propagator 
for xo 7^ 2 / 0 , 


S{p;xo,yo) = 


Tp^ip, xo)Nf{p, 2 / 0)75 for Xo < 2 / 0 , 
Tp^{p, xo)A^'’(p, 2 / 0)75 for xo > 2 / 0 - 


(B.28) 
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The numerical construction of the propagators 


In order to determine and N^, one has to impose the reqnired symmetry of 
the propagator, 

75S'(p;xo, 2 / 0)75 = <S'(p;2/o,xo)^ (B.29) 


which yields 


N^ip,yo) = V(p,2/o)V'^(p,2/o)^ (B.30) 

N\p,yo) = V(p,2/o)V^(p,2/o)^- (B.31) 


Knowing that (|B.18|) is valid especially for zq = xq + 1 and zq = xq — I and 
using (E23). one now concludes that, for the case xq = yo, one has 


S'(p;a;o,a;o) = P-'0^(p, a;o)V(p, xo)V'^(p, a;o)'^75 

+F+'0'’(p, a;o)V(p, Xq) V^(P, (B.32) 


Now, one is left with the task to determine V(p, Xq). A lengthy calculation shows 
that one has 


[V(p,xo)^] ^ = [P_V'^(p,xo)]'^75^+V''’(p,2;o - 1) 

-[P+'iIj^{p,Xq - 1)]H5 P_?/)'’(p,Xo). 


Furthermore, this expression can be shown to be independent of Xq, 


(B.33) 


V(p,xo) = V(p). 


(B.34) 


( |B.33D may hence be evaluated at an arbitrarily chosen xq. Due to the boundary 
conditions, the expression becomes particularly simple for xq = 1, where one has 


[V(p)T ^ = [P_V'-^(p, 1)]H5P+V'^P,0). 


Having computed V, one can now construct the whole propagator. 


>5(p;a^o,2/o) 


Xo)V(p)'0'’(p, 2 / 0 )^5 for Xo < 2 / 0 , 

'^'’(p, Xo)V(p) V^(p, 220)^5 for Xq > yo, 

Xo)V(p)?/^^(p, Xo)^75 

+P+i>\p, Xo)V{p)^i>-^{p, Xo)h5 for Xo = yo- 


(B.35) 


(B.36) 






Appendix C 

The extrapolation procedure 


For the extrapolation of the 1-loop data to the continuum limit, the method 


described in 
used in H75 


5^ was used, which is a generalisation of the blocking technique 


In order to simplify the notation, the following description is given in lattice 
units, i. e. a = 1. The results which are to be extrapolated are n numbers F[L) 
for a given range Li < L 2 ... < Ln- The aim is to extract the leading coefficient 
for L —> 00 . 

One source of errors are roundoff effects. Optimistically, one may assume that 
the values F{L) are correct up to machine precision. 


5p{L) = e\F{L)l (C.l) 

with e ~ 10“^^ for double precision arithmetic, which was used for all compu¬ 
tations in this thesis. Since the calculation of the Feynman diagrams involves 
large sums of terms of different sign, this is clearly an underestimation. In order 
to get a more realistic idea of the size of the roundoff errors, one could compare 
double precision to single (or quadruple) precision results and estimate the L- 
dependence of e. This was done, for example, with the 2-loop results for the 
coupling in [^], where one found a growth of the roundoff error proportional 
to for most contributions. However, experience shows that the error of the 
extrapolation usually is completely dominated by the systematic error, so that 
this procedure for estimating the roundoff error did not seem to be necessary for 
our calculations. 

For the extrapolation, we assume F{L) to have an asymptotic expansion in 
functions fk{L) with k = 1,2,..., Uf, 


nj 


F{L)=J2(^kfk{L)+R{L), 


k=l 


where the rest R{L) behaves like 

R(L) 


k,(L) 


0 as L —> cx). 


(C.2) 


(C.3) 
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The extrapolation procedure 


For improved 1-loop quantities, one would choose the functions fk{L) to be 
l,lnL/L^, l/L^,lnL/L^, 1/L^.... Writing the n data values F{L) as an n- 
dimensional column vector, (|C.2|) becomes 


F = fa + R, 


(C.4) 


where / is an n x n/ matrix and a is the Uf dimensional vector one wants to 
determine. 

In order to get a, one has to minimise the quadratic form 

X^ = {F-fayW\F-fa). (C.5) 

Here, is an n x n matrix of positive weights. It can be used to put an emphasis 
on small or large L. While the values at small L are less affected by roundoff 
errors, the asymptotic expansion is expected to hold to a better degree at large 
L, resulting in a smaller systematic error. In this thesis, all extrapolations where 
done with W = 1, and the behaviour of the systematic errors was examined by 
taking an L range from to F^ax and varying Lmin at fixed F^ax- However, 
for reasons of generality, the procedure will be outlined here for arbitrary W^. 

Minimisation of (|C.5|) yields 

fW^fa = fW^F. (C.6) 


The columns of Wf, which are assumed to be linearly independent, span an 
n/-dimensional subspace. Denoting the projector onto this subspace by P, one 
gets 

Wfa = PWF. (C.7) 


This equation has to be solved for a. As a suitable way to do this, the singular 
value decomposition for Wf turned out to give stable results. To this end, 
one uses the factorisation 

Wf = USV^, (C.8) 


where f/ is an n x u/ matrix with orthonormal columns and 


U^U = 1, UU^ = P. 


(C.9) 


S and V are both Uf x rif matrices, S being diagonal and V orthonormal. Thus 
one gets the solution for a, 

a = VS-^U^WF. (C.IO) 


Now, one has to determine the error of ak- The roundoff error 6a^. is obtained 
from (|C.1| ) by simple error propagation. 


(C.ll) 











Ill 


The determination of the systematic error is a more delicate problem. In 
it was found that a convenient method is the following. 

One assumes that the remainder R can be modelled by a linear combination 
of Ur functions /n^+i, • • •, fnf+rir- the 1-loop calculations in this thesis, one 
may choose the functions 1, InL/L, 1/L,..., InL/L™', 1/L”^ with some positive 
integer m for the fit and take the riy = 2 functions InL/L"*"*"^ and to 

model the remainder. Now, one does separate fits including the Uf functions 
used before and one of the extra functions. In these fits, one gets coefficients 
Ai,... An^ for the extra functions. In order to estimate the error, one repeats 
the original fit, this time not using the data F{L) but the function Aifnf+i{L). 
In this fit, each of the nj functions fk will get a coefficient j3k- The same is done 
with the other — 1 extra functions. For each function f^, one thus gets 
coefficients l3k- The largest of these jSk is then taken as the systematic error 
of the coefficient 
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Appendix D 

Tables of expansion coefficients 


The 1-loop coefficient which is set to zero in order to compute the critical 
bare mass, can be decomposed into 


= ui + NiU2 + + cX^’u^ + c^’u^ + arriX^u^^ + cX’ut. 


di) 


:(!)„ 




,(b. 


di)„ 


(D.l) 


The coefficients Ui are shown in tables p.l| and p.2| for 6 = Ti/b and p = 1 and 
in tables p.3| and p.4| for the other values of 9 and p. 

In the same way, and m' (defined using the upper boundary quark 

fields) may be decomposed according to 

am{‘ (2L/ a) = Vi + NiV 2 + + arrXX^ (L/a)n6 + (D.2) 

with the coefficients Vi for 9 = ir/b and p = 1 in tables and p.6| , and 

=wi + NiW2 + c^w^a + + amX^WQ + (D-3) 

with the coefficients Wi for 9 = and p = 1 in tables D.7 and D.8. 


The renormalisation constant Zp has the coefficients 
2p ’ = Zi + NiZ2 + + cP^4 + 


(DA) 


where the coefficients Zi can be found in tables OT and D.IO for 9 = ii/b and 
p = 1 and in tables p.ll| and p.l2| for different values of 9 and p. Finally, 
the renormalisation constant Zp at 2L such that mi{L) = 0 may be expanded 
according to 

Z^^^(2L/a) = zi + N{Z2 + 4i^Z3 + c^^^za + amX^(L/a)z6 + 4^^Z7, (D.5) 

with the coefficients Zi in tables p.l3| and p.l4| for 9 = 7r/5 and p = 1 and in 
tables p.l5| and p.l6| for the other values of 9 and p. 

In the tables, up to ten digits of the coefficients are shown. Since they have 
been computed using double precision arithmetic and experience shows that up 
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Tables of expansion coefficients 


to L/a = 32, one looses two to three digits due to roundoff errors, all digits shown 
here should be significant. To get a better estimate for the roundoff errors, one 
would have to compare the results for some small lattices with single or quadruple 
precision results and extract the scaling behaviour of the errors for increasing 
L/a. In this thesis, however, we are only interested in discretisation errors, where 
the last digits are irrelevant, and in extrapolations, which are dominated by the 
systematic errors due to higher order terms in the expansion. Therefore, a precise 
determination of the roundoff errors did not seem to be necessary. 


L/a 

Ml 

U2 

U3 

4 

0.3602839939 

-0.003225474813 

-0.09605428870 

5 

0.3232002646 

-0.000774641219 

-0.05739586792 

6 

0.3037717960 

-0.000214828175 

-0.03886722333 

7 

0.2945227639 

-0.000090081829 

-0.02788612762 

8 

0.2886435020 

-0.000047062385 

-0.02116169433 

9 

0.2846484024 

-0.000027131359 

-0.01653724923 

10 

0.2818166195 

-0.000017284660 

-0.01334401178 

11 

0.2797403581 

-0.000011720940 

-0.01096120300 

12 

0.2781681231 

-0.000008284863 

-0.00919331010 

13 

0.2769535573 

-0.000006129109 

-0.00780391809 

14 

0.2759917826 

-0.000004608993 

-0.00672232221 

15 

0.2752204783 

-0.000003596909 

-0.00584115805 

16 

0.2745896768 

-0.000002820717 

-0.00513109826 

17 

0.2740694316 

-0.000002283224 

-0.00453708978 

18 

0.2736334509 

-0.000001846627 

-0.00404578595 

19 

0.2732659914 

-0.000001535330 

-0.00362632773 

20 

0.2729521158 

-0.000001271725 

-0.00327224614 

21 

0.2726829531 

-0.000001079336 

-0.00296502249 

22 

0.2724494848 

-0.000000911236 

-0.00270139289 

23 

0.2722464289 

-0.000000786142 

-0.00246962788 

24 

0.2720680657 

-0.000000674154 

-0.00226804856 

25 

0.2719111055 

-0.000000589418 

-0.00208889322 

26 

0.2717717705 

-0.000000512088 

-0.00193130204 

27 

0.2716479341 

-0.000000452721 

-0.00178994940 

28 

0.2715370145 

-0.000000397695 

-0.00166441225 

29 

0.2714375933 

-0.000000354905 

-0.00155092337 

30 

0.2713478540 

-0.000000314740 

-0.00144929720 

31 

0.2712668249 

-0.000000283144 

-0.00135679872 

32 

0.2711931957 

-0.000000253175 

-0.00127337243 


Table D.l: Parts of the 1-loop coefficient for 6 = n/5 and p = 1 
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L/a 

U4 

Ms 

Uq 

Uj 

4 

0.7316730490 

0.6084363759 

1.188647276 

0.01803922007 

5 

0.3210591192 

0.3878385759 

1.117262157 

0.00413349286 

6 

0.0042189812 

0.2631884549 

1.079778023 

-0.00157603839 

7 

0.0024050333 

0.1935221158 

1.057964363 

-0.00069723702 

8 

0.0012919050 

0.1464623871 

1.043886577 

-0.00038756353 

9 

0.0008286417 

0.1158988293 

1.034483484 

-0.00020611832 

10 

0.0005124163 

0.0932417601 

1.027741909 

-0.00012753286 

11 

0.0003555287 

0.0771715444 

1.022853029 

-0.00007616582 

12 

0.0002404411 

0.0645591187 

1.019116087 

-0.00005074611 

13 

0.0001763427 

0.0550783445 

1.016254342 

-0.00003280704 

14 

0.0001268701 

0.0473445009 

1.013970372 

-0.00002306893 

15 

0.0000969510 

0.0412863218 

1.012152254 

-0.00001579511 

16 

0.0000729793 

0.0362045275 

1.010655337 

-0.00001157037 

17 

0.0000575581 

0.0320993397 

1.009428779 

-0.00000827330 

18 

0.0000448460 

0.0285822227 

1.008394882 

-0.00000625754 

19 

0.0000362650 

0.0256723716 

1.007528515 

-0.00000462840 

20 

0.0000290320 

0.0231377112 

1.006784665 

-0.00000359173 

21 

0.0000239575 

0.0210003357 

1.006150100 

-0.00000272952 

22 

0.0000196031 

0.0191135425 

1.005597146 

-0.00000216309 

23 

0.0000164500 

0.0174975213 

1.005118494 

-0.00000168054 

24 

0.0000137045 

0.0160552118 

1.004696293 

-0.00000135523 

25 

0.0000116632 

0.0148037736 

1.004326361 

-0.00000107234 

26 

0.0000098641 

0.0136765405 

1.003996718 

-0.00000087756 

27 

0.0000084961 

0.0126877084 

1.003704907 

-0.00000070513 

28 

0.0000072780 

0.0117900311 

1.003442618 

-0.00000058430 

29 

0.0000063338 

0.0109951461 

1.003208384 

-0.00000047564 

30 

0.0000054857 

0.0102686695 

1.002996274 

-0.00000039839 

31 

0.0000048172 

0.0096201313 

1.002805407 

-0.00000032793 

32 

0.0000042121 

0.0090239315 

1.002631448 

-0.00000027722 


Table D.2: Parts of the 1-loop coefficient for 6 = Ti/h and p = 1 (cont.) 
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L/a 

Ml 

U2 

U3 

9 = 0, p=l 

4 

0.3117291700 

-0.001413989874 

-0.06325439280 

6 

0.2851814233 

-0.000072348720 

-0.02653692904 

8 

0.2786313238 

-0.000015427643 

-0.01462634537 

10 

0.2755771628 

-0.000005279182 

-0.00927252482 

12 

0.2739097717 

-0.000002238838 

-0.00640604311 

14 

0.2729004050 

-0.000001074282 

-0.00469179577 

16 

0.2722432007 

-0.000000559106 

-0.00358486772 

9 = 0.5, p= 1 

4 

0.3527716036 

-0.003134326041 

-0.08882286240 

6 

0.3009793922 

-0.000203934067 

-0.03600623562 

8 

0.2872085538 

-0.000044256296 

-0.01960752669 

10 

0.2809505241 

-0.000016156714 

-0.01236285801 

12 

0.2775900272 

-0.000007687118 

-0.00851637431 

14 

0.2755788098 

-0.000004247030 

-0.00622673813 

16 

0.2742799815 

-0.000002584214 

-0.00475246807 

9 = 0, p = 2 

4 

0.2890122981 

-0.000466276200 

-0.02966169499 

6 

0.2776917924 

-0.000096885717 

-0.01285486035 

8 

0.2742017998 

-0.000032373848 

-0.00717204102 

10 

0.2726668499 

-0.000014117091 

-0.00457353269 

12 

0.2718544505 

-0.000007249144 

-0.00317000300 

9 = 0.5, p = 2 

4 

0.3230287075 

0.004646415684 

-0.04091643990 

6 

0.2928573371 

0.000682358234 

-0.01682956807 

8 

0.2827493514 

0.000177104745 

-0.00920123428 

10 

0.2781432425 

0.000066678908 

-0.00581015860 

12 

0.2756603756 

0.000030966376 

-0.00400514977 


Table D.3: Parts of the 1-loop coefficient rni'^ 


for various values of 6 and p 
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L/a 

U4 

Ms 

Me 

u^ 



0 = 0, p 

= 1 


4 

0.5887651776 

0.2773460352 

1.060808688 

0.01414891282 

6 

0.0007741031 

0.1200813871 

1.027928002 

-0.00039927726 

8 

0.0001985526 

0.0670949370 

1.016024884 

-0.00007249879 

10 

0.0000659113 

0.0428310908 

1.010376100 

-0.00001642414 

12 

0.0000255381 

0.0297073108 

1.007259112 

-0.00000361335 

14 

0.0000108268 

0.0218109132 

1.005360129 

-0.00000027458 

16 

0.0000047721 

0.0166920181 

1.004118703 

0.00000058142 



9 = 0.5, p 

= 1 


4 

0.7136493082 

0.5749737393 

1.168810060 

0.01714141379 

6 

0.0030661542 

0.2495987793 

1.072695582 

-0.00138698280 

8 

0.0009267006 

0.1391978462 

1.040386927 

-0.00033684366 

10 

0.0003656376 

0.0887236823 

1.025677627 

-0.00010975186 

12 

0.0001711788 

0.0614750484 

1.017759437 

-0.00004324006 

14 

0.0000902366 

0.0451034411 

1.013012168 

-0.00001945044 

16 

0.0000518894 

0.0345013971 

1.009943004 

-0.00000964499 



o' 

= 2 


4 

0.0032392824 

0.0936449462 

1.036188030 

-0.00040229454 

6 

0.0005702345 

0.0394953561 

1.014454343 

-0.00004485954 

8 

0.0001630428 

0.0217769672 

1.007693781 

-0.00000952668 

10 

0.0000612800 

0.0138044243 

1.004761213 

-0.00000285504 

12 

0.0000273938 

0.0095357430 

1.003232426 

-0.00000105727 



9 = 0.5, p 

= 2 


4 

0.0011902517 

0.4959636094 

1.137886597 

-0.00272630910 

6 

0.0001742293 

0.2178875276 

1.060855811 

-0.00034621022 

8 

0.0000433483 

0.1220946720 

1.034142311 

-0.00008143737 

10 

0.0000147043 

0.0780067469 

1.021817535 

-0.00002664492 

12 

0.0000060319 

0.0541217346 

1.015134789 

-0.00001071394 


Table D.4: Parts of the 1-loop coefficient for various values of 6 and p (cont.) 
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Tables of expansion coefficients 


L/a 

Vl 

V2 

^3 

4 

0.2888594021 

-0.00005058908107 

-0.02117043086 

5 

0.2820887015 

-0.00001972144488 

-0.01335256814 

6 

0.2782616600 

-0.00000882525879 

-0.00919576341 

7 

0.2760809582 

-0.00000497218858 

-0.00672431079 

8 

0.2746373747 

-0.00000296561503 

-0.00513202789 

9 

0.2736756316 

-0.00000194609439 

-0.00404651418 

10 

0.2729796791 

-0.00000132369711 

-0.00327267396 

11 

0.2724733713 

-0.00000094810144 

-0.00270172931 

12 

0.2720853970 

-0.00000069646118 

-0.00226827213 

13 

0.2717867718 

-0.00000052844286 

-0.00193148046 

14 

0.2715486018 

-0.00000040853487 

-0.00166454014 

15 

0.2713579394 

-0.00000032292580 

-0.00144940101 

16 

0.2712013150 

-0.00000025894756 

-0.00127345074 


Table D.5: Parts of the 1-loop coefficient m^i\2L) with mi(L) = 0 for 6 = vr/S 
and p = 1 


L/a 

Va 

V5 

Vq 

V7 

4 

0.001291461837 

0.1463868229 

1.045074777 

-0.0003869872089 

5 

0.000511742035 

0.0931663422 

1.029247500 

-0.0001269718170 

6 

0.000240252491 

0.0645364689 

1.019635452 

-0.0000506118783 

7 

0.000126734771 

0.0473260174 

1.014466186 

-0.0000229801784 

8 

0.000072924806 

0.0361957674 

1.010920783 

-0.0000115365918 

9 

0.000044809374 

0.0285753372 

1.008629730 

-0.0000062357257 

10 

0.000029013497 

0.0231336427 

1.006938179 

-0.0000035810441 

11 

0.000019590549 

0.0191103368 

1.005730206 

-0.0000021559966 

12 

0.000013697243 

0.0160530751 

1.004792851 

-0.0000013512045 

13 

0.000009858993 

0.0136748331 

1.004080301 

-0.0000008747879 

14 

0.000007274769 

0.0117888051 

1.003507183 

-0.0000005825718 

15 

0.000005483328 

0.0102676735 

1.003052472 

-0.0000003971558 

16 

0.000004210539 

0.0090231792 

1.002676692 

-0.0000002763990 


Table D.6: Parts of the 1-loop coefficient rri{\2L) with mi{L) = 0 for 9 = vr/S 
and p = 1 (cont.) 
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L/a 

Wi 

W2 

W3 

4 

0.2971376852 

0.007159453631 

0.07069378144 

5 

0.2930551099 

0.002002507543 

0.04435489541 

6 

0.2845876866 

0.000904887361 

0.03097531507 

7 

0.2820563435 

0.000384152314 

0.02227017749 

8 

0.2796089169 

0.000228377566 

0.01706219846 

9 

0.2779500101 

0.000127950492 

0.01329286954 

10 

0.2764995647 

0.000087862653 

0.01077351137 

11 

0.2755037279 

0.000057206374 

0.00882189295 

12 

0.2746365666 

0.000042530198 

0.00741823947 

13 

0.2740124775 

0.000030217871 

0.00628082481 

14 

0.2734635897 

0.000023601409 

0.00541982387 

15 

0.2730515248 

0.000017755117 

0.00469970451 

16 

0.2726849153 

0.000014352239 

0.00413374337 

17 

0.2724001308 

0.000011245712 

0.00364919780 

18 

0.2721440627 

0.000009327024 

0.00325731535 

19 

0.2719396253 

0.000007535561 

0.00291573391 

20 

0.2717540932 

0.000006376663 

0.00263317208 

21 

0.2716026412 

0.000005276671 

0.00238335648 

22 

0.2714640928 

0.000004537982 

0.00217290423 

23 

0.2713489030 

0.000003828045 

0.00198468504 

24 

0.2712427988 

0.000003336311 

0.00182372283 

25 

0.2711532200 

0.000002859114 

0.00167838569 

26 

0.2710702125 

0.000002519843 

0.00155251931 

27 

0.2709992153 

0.000002188062 

0.00143795523 

28 

0.2709330829 

0.000001946840 

0.00133766988 

29 

0.2708758835 

0.000001709481 

0.00124576235 

30 

0.2708223594 

0.000001533525 

0.00116456364 

31 

0.2707756142 

0.000001359513 

0.00108970680 

32 

0.2707316956 

0.000001228300 

0.00102303872 


Table D.7: Parts of m' at 1-loop order with mi(L) = 0 for 6 = tt/ 5 and p = 1 
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Tables of expansion coefficients 


L/a 

W4 

W5 

Wq 

Wj 

4 

0.6670524594 

0.8104561493 

1.111943322 

-0.01095404633 

5 

0.3903966722 

0.5388856057 

1.091349905 

-0.00357819396 

6 

-0.0022730281 

0.3747386559 

1.069419226 

0.00063957047 

7 

-0.0010901031 

0.2781810491 

1.055379436 

0.00017422876 

8 

-0.0006270025 

0.2119249521 

1.043698522 

-0.00008757089 

9 

-0.0003743296 

0.1681735583 

1.035783207 

-0.00008690205 

10 

-0.0002493311 

0.1356020041 

1.029322006 

-0.00009872911 

11 

-0.0001683210 

0.1123440306 

1.024688617 

-0.00007582581 

12 

-0.0001217515 

0.0940684807 

1.020845628 

-0.00006576171 

13 

-0.0000885627 

0.0802868921 

1.017957986 

-0.00005052598 

14 

-0.0000674800 

0.0690416443 

1.015515375 

-0.00004195232 

15 

-0.0000515736 

0.0602178847 

1.013611180 

-0.00003293300 

16 

-0.0000407400 

0.0528164319 

1.011971652 

-0.00002726629 

17 

-0.0000322524 

0.0468314608 

1.010655815 

-0.00002188842 

18 

-0.0000261695 

0.0417044433 

1.009505714 

-0.00001826479 

19 

-0.0000212731 

0.0374600536 

1.008560908 

-0.00001495738 

20 

-0.0000176252 

0.0337634280 

1.007724628 

-0.00001261428 

21 

-0.0000146275 

0.0306449645 

1.007024441 

-0.00001050889 

22 

-0.0000123250 

0.0278924569 

1.006398071 

-0.00000895924 

23 

-0.0000104015 

0.0255343162 

1.005865307 

-0.00000757449 

24 

-0.0000088871 

0.0234298978 

1.005384399 

-0.00000652389 

25 

-0.0000076047 

0.0216036230 

1.004969907 

-0.00000558578 

26 

-0.0000065743 

0.0199587646 

1.004592861 

-0.00000485624 

27 

-0.0000056916 

0.0185156623 

1.004264203 

-0.00000420373 

28 

-0.0000049703 

0.0172056983 

1.003963238 

-0.00000368576 

29 

-0.0000043461 

0.0160456199 

1.003698336 

-0.00000322116 

30 

-0.0000038285 

0.0149854514 

1.003454342 

-0.00000284588 

31 

-0.0000033768 

0.0140389494 

1.003237764 

-0.00000250813 

32 

-0.0000029975 

0.0131688820 

1.003037259 

-0.00000223120 


Table D.8: Parts of m' at 1-loop order with mi(L) = 0 for 6 = tt/ 5 and p = 1 
(cont.) 
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L/a 


^2 

Z3 

4 

0.09550891946 

0.03423528831 

-0.1341134935 

6 

0.16565857330 

0.01473258948 

-0.0932504829 

8 

0.22869553805 

0.00949691345 

-0.0711779334 

10 

0.28916821812 

0.00753246481 

-0.0574507258 

12 

0.34887762674 

0.00654228932 

-0.0481191321 

14 

0.40849882024 

0.00594073122 

-0.0413752997 

16 

0.46828819102 

0.00553322566 

-0.0362791400 

18 

0.52833626903 

0.00523794769 

-0.0322950506 

20 

0.58866515448 

0.00501397159 

-0.0290961208 

22 

0.64926855814 

0.00483827767 

-0.0264717741 

24 

0.71012902785 

0.00469683169 

-0.0242803679 

26 

0.77122555828 

0.00458056303 

-0.0224231782 

28 

0.83253695525 

0.00448334105 

-0.0208293183 

30 

0.89404326241 

0.00440087155 

-0.0194465900 

32 

0.95572628769 

0.00433005700 

-0.0182357146 


Table D.9: Parts of the 1-loop coefficient Zp^ 


at 6* = 7r/5 and p 


1 


L/a 

^4 

^6 

Z7 

4 

0.2752077673 

1.240156503 

-0.1954547977 

6 

0.3024134295 

1.582932389 

-0.1615028450 

8 

0.2903043621 

1.870599424 

-0.1331757885 

10 

0.2677764380 

2.136141388 

-0.1121302285 

12 

0.2446013168 

2.391341118 

-0.0964095495 

14 

0.2234839233 

2.641100707 

-0.0843763356 

16 

0.2049375164 

2.887730110 

-0.0749275370 

18 

0.1888173394 

3.132429195 

-0.0673363702 

20 

0.1748109198 

3.375870929 

-0.0611160764 

22 

0.1625955840 

3.618456709 

-0.0559323228 

24 

0.1518850825 

3.860438271 

-0.0515494198 

26 

0.1424386316 

4.101980209 

-0.0477971213 

28 

0.1340577488 

4.343193976 

-0.0445497094 

30 

0.1265798562 

4.584157332 

-0.0417124957 

32 

0.1198717596 

4.824925943 

-0.0392128935 


Table D.IO: Parts of the 1-loop coefficient Zp^ at 0 = tt/S and p = 1 (cont.) 
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Tables of expansion coefficients 


L/a 


^2 

^3 

9 = 0, p=l 

4 

0.1326808072 

0.01528384278 

-0.1168458223 

6 

0.1474652524 

0.00728806063 

-0.0779791840 

8 

0.1731329290 

0.00439645001 

-0.0585101791 

10 

0.2035733306 

0.00321873110 

-0.0468218991 

12 

0.2367057252 

0.00262320959 

-0.0390261029 

14 

0.2715872938 

0.00226224002 

-0.0334556609 

16 

0.3077062797 

0.00201667075 

-0.0292766717 


9 

= 0.5, p = 1 


4 

0.0824520269 

0.03143610800 

-0.1238972498 

6 

0.1271965066 

0.01394010853 

-0.0857279226 

8 

0.1683678316 

0.00887347640 

-0.0652575226 

10 

0.2089123772 

0.00693665488 

-0.0525900075 

12 

0.2497370083 

0.00597030160 

-0.0440062052 

14 

0.2910727552 

0.00538995423 

-0.0378156722 

16 

0.3329441002 

0.00499976060 

-0.0331442026 

9 = 0, p = 2 

4 

0.6521364999 

0.02993994420 

-0.1209100243 

6 

0.9221596868 

0.01125513464 

-0.0810128062 

8 

1.1964703064 

0.00632783134 

-0.0608921086 

10 

1.4728928412 

0.00424833271 

-0.0487675434 

12 

1.7507932870 

0.00311330602 

-0.0406652864 


0 

= 0.5, p = 2 


4 

0.0665300894 

0.16011586673 

-0.0811127180 

6 

0.1791766558 

0.07550278985 

-0.0548862358 

8 

0.2766417147 

0.04877505559 

-0.0413683000 

10 

0.3693268329 

0.03915526808 

-0.0331678255 

12 

0.4605318775 

0.03486305919 

-0.0276723609 


Table D.ll: Parts of the 1-loop coefficient Zp^ 


at various values of 6 and p 





























123 


L/a 

Z4 

^6 

Z7 


9 

= 0, P = 1 


4 

0.2448721196 

1.297122470 

-0.1526509731 

6 

0.1942219137 

1.476291125 

-0.1197969999 

8 

0.1640794499 

1.643098574 

-0.0965623402 

10 

0.1420781526 

1.806546577 

-0.0803607200 

12 

0.1251135646 

1.968711620 

-0.0686367894 

14 

0.1116498521 

2.130300176 

-0.0598239325 

16 

0.1007298970 

2.291603872 

-0.0529813561 


9 -- 

= 0.5, p= 1 


4 

0.2334171022 

1.189005361 

-0.1834591369 

6 

0.2443514694 

1.448427447 

-0.1510297236 

8 

0.2314225501 

1.659729841 

-0.1241302945 

10 

0.2123099369 

1.853796129 

-0.1042991779 

12 

0.1934250942 

2.040359892 

-0.0895616671 

14 

0.1764740188 

2.223182810 

-0.0783182317 

16 

0.1616952383 

2.403947754 

-0.0695088918 


9 

= 0, p = 2 


4 

0.5124378929 

2.975968192 

-0.1883730541 

6 

0.4114552065 

4.088045527 

-0.1353705625 

8 

0.3378163059 

5.172972867 

-0.1049844154 

10 

0.2849616955 

6.247086041 

-0.0855561121 

12 

0.2458461197 

7.315968086 

-0.0721305763 


9 = 

= 0.5, p = 2 


4 

0.2819915539 

1.490852756 

-0.2128630913 

6 

0.2895835015 

1.938807978 

-0.1612690338 

8 

0.2645149383 

2.334410619 

-0.1275581165 

10 

0.2365500613 

2.709957987 

-0.1049282241 

12 

0.2116794602 

3.076362839 

-0.0889178266 


Table D.12: Parts of the 1-loop coefficient Zp^ at various values of 9 and p (cont.) 
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Tables of expansion coefficients 


L/a 


^2 


4 

0.2293123742 

0.009573101967 

-0.07147065652 

6 

0.3494732430 

0.006579092370 

-0.04825027064 

8 

0.4688244330 

0.005555581838 

-0.03634661389 

10 

0.5891478869 

0.005029207837 

-0.02913523296 

12 

0.7105656939 

0.004707933201 

-0.02430499001 

14 

0.8329343379 

0.004491803200 

-0.02084578712 

16 

0.9560901651 

0.004336724846 

-0.01824725568 


Table D.13: Parts of the 1-loop coefficient Zp\2L) with mi{L) = 0 at 6 = tt/ 5 
and p = I 


L/a 

^4 

^6 


4 

0.2911905373 

1.874547330 

-0.1335040402 

6 

0.2450740029 

2.394596702 

-0.0965677693 

8 

0.2051995671 

2.890423111 

-0.0750119316 

10 

0.1749691162 

3.378169482 

-0.0611660627 

12 

0.1519872262 

3.862442505 

-0.0515813433 

14 

0.1341272577 

4.344969442 

-0.0445712831 

16 

0.1199210816 

4.826518506 

-0.0392281295 


Table D.14: Parts of the 1-loop coefficient Zp\2L) with mi{L) = 0 at 9 = n/h 
and p = I (cont.) 
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L/a 


^2 

h 

9 = 0, p=l 

4 

0.1844808442 

0.004372112061 

-0.05984829671 

6 

0.2453504776 

0.002601676822 

-0.03949477059 

8 

0.3142628791 

0.002001926120 

-0.02947930729 

10 

0.3877224203 

0.001689952794 

-0.02352815290 

12 

0.4639602752 

0.001496341488 

-0.01958158454 

14 

0.5420480669 

0.001364127036 

-0.01677124455 

16 

0.6214500817 

0.001268049176 

-0.01466752787 

9 = 0.5, p = 1 

4 

0.1689995368 

0.008923199948 

-0.06548732132 

6 

0.2505009004 

0.005999143948 

-0.04413591317 

8 

0.3336383383 

0.005016999229 

-0.03321155407 

10 

0.4187398577 

0.004517026013 

-0.02660691792 

12 

0.5054112010 

0.004213452330 

-0.02218859089 

14 

0.5932918622 

0.004009910129 

-0.01902679794 

16 

0.6821137770 

0.003864208278 

-0.01665284161 

9 = 0, p = 2 

4 

1.2060785837 

0.007491917291 

-0.06149659905 

6 

1.7559778103 

0.003528122640 

-0.04081983117 

8 

2.3131197146 

0.002140383509 

-0.03057928797 

10 

2.8738200566 

0.001463496442 

-0.02445228153 

12 

3.4366941330 

0.001073113077 

-0.02037239719 

9 = 0.5, p = 2 

4 

0.2773639094 

0.049707671735 

-0.04185048543 

6 

0.4609370909 

0.035128446020 

-0.02778302053 

8 

0.6423873199 

0.031226685817 

-0.02081561802 

10 

0.8244441384 

0.029527745468 

-0.01664726569 

12 

1.0074272069 

0.028606927366 

-0.01387127057 


Table D.15: Parts of the 1-loop coefficient Zp\2L) with mi{L) = 0 at various 
values of 9 and p 
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Tables of expansion coefficients 


L/a 

^4 

^6 


9 = 0, p=l 

4 

0.1718305672 

1.687871282 

-0.09775258579 

6 

0.1280139314 

2.002040662 

-0.06907686487 

8 

0.1020232090 

2.316619630 

-0.05317715374 

10 

0.0848215601 

2.633590801 

-0.04317899772 

12 

0.0725826258 

2.952117432 

-0.03632993149 

14 

0.0634269015 

3.271619233 

-0.03134989520 

16 

0.0563194892 

3.591757515 

-0.02756760705 



= 0.5, p= 1 


4 

0.2321338155 

1.663348371 

-0.12438362939 

6 

0.1938923565 

2.044179271 

-0.08971537396 

8 

0.1619544256 

2.407180440 

-0.06959163712 

10 

0.1379679076 

2.765026025 

-0.05670609836 

12 

0.1198008298 

3.120873229 

-0.04780137447 

14 

0.1057072432 

3.475803721 

-0.04129513010 

16 

0.0945069004 

3.830264175 

-0.03633907736 

9 = 0, p = 2 

4 

0.3405246084 

5.207874762 

-0.10562217877 

6 

0.2465730856 

7.335024524 

-0.07229922856 

8 

0.1926801273 

9.459097330 

-0.05491154830 

10 

0.1579528577 

11.581344073 

-0.04425113515 

12 

0.1337701693 

13.702563787 

-0.03705152726 


9 

= 0.5, p = 2 


4 

0.2661600367 

2.343188522 

-0.12821614899 

6 

0.2121269758 

3.079864030 

-0.08907791880 

8 

0.1730920116 

3.798786095 

-0.06801232921 

10 

0.1454143632 

4.511560279 

-0.05494446490 

12 

0.1251028064 

5.221609660 

-0.04606774596 


Table D.16: Parts of the 1-loop coefficient Zp\2L) with mi{L) = 0 at various 
values of 9 and p (cont.) 
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